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NATO Advanced Study Institute: 
Granular Nanoelectronics 


July 23.— August 4, 1990, i Ciocco, Italy 


Please remember that the role of the Advanced-Study Institute is one of education. The talks 
reflect some introductory material and make connection with the.other speakers, while at the 
same time presenting up-to-date research. Each-will try to relate in some way to the overall 
theme of the ASI, which is granular nanoelectronics — the future of few-electron systems for 
electronics applications. Information is best transmitted by a lively question and answer 
period. We encourage each attendee to ask questions and make sure that the topic is made 
understandable to him (or her). 


Talks 1-3 are morning talks each day, while talks 4,5 are afternoon talks. 


Program 
Monday, 23 July 
Opening Remarks and general welcome 


quantum systems" 

Dave Ferry, "Projections on the future of ULSI and nanolithograpy” 
Steve Beaumont, “Nanofabrication of quantum wires and rings" 
Karl Hess, "Path integrals and transport in mesoscopic systems” 


Tuesday, 24 July 


ta BR WD Nr 


Markus Biittiker, "Transmission probabilities and electric resistance" 

Alan Fowler,"Interference phenomena in devices structures" 

John-Barker, "Introduction to quantum transport in quantum waveguides" 
Tony Leggett, “Dephasing and non-dephasing collisions in nanostructures" 
Trevor Thornton, "Electronic fluctuations and transport in waveguides" 


Wednesday, 25 July 


ARON 


Jérg Kotthaus, "Spectroscopy of electronic excitations in low dimensional systems" 
Alan Fowler, "Waveguide devices: rings, non-local transport, etc." 

POSTER SESSION 1 

Mark Reed, "Tunneling investigations of low dimensional systems" 

Erick Gomik, "Tunneling between constrained dimensionality systems" 


Thursday, 26 July 
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1. Jérg Kotthaus, “Lateral confinement in surface superlattices" 
2. Dave Ferry, "Magnetoconductance in surface superlattices" 
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Carlo Beenakker, "Ballistic and adiabatic transport 1" 
Markus Biittiker, “Quantum oscillations in loops and tunnel junctions" 
Steve Beaumont, "Optical properties of arrays of microdevices” 


riday, 27 July 


Carlo Jacoboni, "Numerical approaches to non-equilibrium quantum transport" 
Mark Reed, "Non-equilibrium quantum effects — structures and devices” 

Lino Reggiani, “Noise in small and ultra-small device structures" 

Gerhardt Abstreiter, “Growth engineering of mesoscopic systems” 

Antti Jauho, "Non-equilibrium Green's.functions in quantum transport" 


\ Monday, 30 July 


Lex Akers, "Neural and constrained interconnect automata" 

J. P. Launay, “Molecular electronics" 

John Barker, "Granular nanoelectronic phenomena" 

Gerhardt Abstreiter, "Quantum field effects on optical properties of superlattices" 
Chihiro Hamaguchi, "Optical properties of short period superlattices” 


Tuesday, 31 July 


Trevor Thornton, "Resistance fluctuations and quantization in microdevices" 
Aaron Szafer, "Fluctuations and quantum chaos in ballistic transport" 

Carlo Beenakker, “Ballistic and adiabatic transport I" 

Friedl Kuchar, "Microwave studies of quasi-one dimensional wires" 
Chihiro Hamaguchi, "Non-equilibrium carrier transport in small structures” 


Wednesday, 1 August 


Aaron Szafer, "Quantum circuits and non-locality" 

Tony Leggett, "Some considerations related to the quantization of charge in mesoscopic 
systems" 

POSTER SESSION 2 

K. K. Likharev, "Single-electronics: The correlated transfer of single electrons in 
ultrasmall tunnel junctions, arrays, and systems" 

L. J. Geerligs, "Charging effects and ‘turnstile’ clocking of single electrons in small 
tunnel junctions” 


Thursday, 2 August 


Toshiaki Ikoma, "Weak localization and phase-breaking mechanisms of electron waves 

in quasi-one-dimensional wires" 

Karl Hess, “Numerical simulations of electron waveguides" 

Lex Akers, “VLSI implementations of neural networks" 

Antti Jauko, "High field quantum transport" 
Lino Reggiani, "Monte Carlo algorithms for non-equilibrium quantum transport" i 


Friday, 3 August 


1. 
2 


Gerry Iafrate, “Correlation and exchange in single electron events" 
Bob O'Connell, "The few-body problem in nanoelectronics" 


John Barker, “Interfacing molecular and biological systems with electronic mesoscopic 
systems" 


CONFERENCE WRAPUP 


Ef 


Poster Manuscripis 


Poster presentations are to last for 3-4 minutes, and should be viewed as an "advertisement" 
for the poster itself. In general, no questions will be entertained during the poster "shotgun" 
presentation. A decision will be made by the organizers during the ASI as to which posters 
should be included in the published proceedings. 


Posters that are accepted for the proceedings will be allowed 4 single-spaced A4 pages (typing 
space of 6.25"x10,194"), which includes all references and figures (manuscripts which 
exceed this limit will be arbitrarily shortened by the editors). The format can be taken from 
the lecture notes provided to attendees, for which the text is prepared originally as 12 point, 
Times-Roman. Original line drawings can be submitted with the manuscript, and will be 
scanned into the document itself. The manuscript is due to’Prof. D. K. Ferry, Department of 
Electrical Engineering, Arizona State University, Tempe, AZ 85287-5706, by September 1, 
1990. A 3.5" disk containing the manuscript as a word processing .ile should be submitted 
with the manuscript. Either Macintosh or MSDOS format is acceptable, although the former 
is preferred. Please note that the full titles and author lists are required for the cited 
references. Equations should be set with MacEgn (preferred to Expressionist), or written in 
by hand for resetting by the editors. 
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Abstracts 


Session 1 
Wednesday, July 25 


1-1 25 to 85nm Gate Length GaAs 
MESFETs and HEMTs, J. M. Ryan, J. Han, and 
D. K. Ferry, Arizona State University, U.S.A. GaAs 
MESFETS and HEMTs with gate lengths from 85 to 
25nm have been fabricated using electron beam 
lithography. Maximum transconductances of 108 
mS/mm and 215 mS/mm were obtained for the 
MESFETs and HEMTs, respectively. These 
apparently low values are in the expected range for 
devices with gate-length-to-channel ratios much 
smaller than oné, so that capacitance fringing is 
important. The transconductance demonstrates 
velocity overshoot is for devices with gate lengths 
below ~55nm. For gates shorter than ~35nm, the 
overshoot-related increase in transconductance 
saturates, This effect is explainable in terms of 
carrier heating in the ohmic contact, which decreases 
the effect of overshoot at the gate. 

The high-frequency response of the devices was 
tested to 20 GHz. The MESFETs yielded a 
maximum cutoff frequency f, of 167 GHz, with f; 
following a linear relation with the inverse gate 
length for gate lengths down to 55 nm. The HEMTs 
had a maximum f, of 18 GHz, a number limited by 
the extra gate capacitance needed for the measurement 
configuration. 

For HEMTs with gate lengths shorter than 30nm, 
the current is found to vary exponentially with the 
channel voltage, Simple calculations indicate that 
this would result when tunneling is the dominant 
current mechanism, Similar behavior has been 
observed in short-gate devices of other groups. For 
gate lengths larger than 30nm, we find the usual 
voltage dependence, For MESFETs with gate lengths 
shorter than 30nm, no evidence of an exponential 
relationship is found, This is probably due to weaker 
confinement of carriers in the MESFETs as compared 
to the HEMTs. 


1-2 Quantitative Determinationof the 
Magnetic-Field Density of States as a 
Function of Fermi Energy in the Two- 
Dimensional Electron Gas, Raymond Ashoori, 
Cornell University, U.S.A. Using capacitive 
techniques, we have studied both the tunneling! 
density of states (DOS) and the thermodynamic DOS 
in a two-dimensional (2D) electron gas of variable 
density (0 to 6x10!!cm=-2) in the presence of 
magnetic field applied perpendicular to the plane of 
the electron gas. Our sample design allows for 
calculation of the thermodynamic DOS as a function 
of Fermi energy of electrons in the 2D gas from 


capacitive data using no smaple paraméters other than 
the sample area. This type of measurement is distinct 
from other DOS measurements which have typically 
measured the DOS as a function of magnetic field at 
fixed density. The DOS determination presented here 
allows for direct comparison with models of the DOS 
as a function of energy. In contrast with fixed-density 
DOS results? which are fit well by Landau levels 
taken to be Gaussian (asa functionof electronic 
energy) plus a background DOS, our results for the 
DOS as a function of Fermi energy are fit best (at 4T 
and lower) with Lorentzian lineshapes. Further, at 
high fields (8.5T) the exchange-cnhanced spin 
splitting is observed, and the exchange energy for 
electrons is determined. Finally, methods and results 


of application of these capacitive techniques to | 


quantum dot arrays patterned on our samples will also 
be presented briefly. 

1 R. C. Ashoori, J. A. Lebens, N. P. Bigelow, and 
R. H. Silsbee, Phys. Rev. Lett. 64, 681 (1990). 

2 See for example E. Gornik, R. Lassnig, G. 
Suasser, H. L. Stérmer, A. C. Gossard, and W. 
Weigmann, Phys. Rey. Lett. 54, 1820 (1985). 


1-3 Study of Lower-Dimensional Transport 
by Electroluminescence, Hans P. Zappe, 
Fraunhofer Institut fiir Angewandte F estkérperphysik, 
West Germany. It is well known that many 
semiconducting materials and devices emit visible 
light when subject to high electric fields, This 
electroluminescence is becoming apparent at quite 
low operating voltages in modern, small-sized 
structures and may be used to study the transport 
behavior in transistors as well as 2- and 1- 
dimensional conducting layers. We have used this 
effect to study the scattering, real-space transfer and 
energy-loss mechanisms of conduction electrons in 
some of these systems, 

The voltage-dependent intensity, lateral distribution 
and energy spectrum of the electroluminescence has 
been measured in a variety of structures. Thee nature 
of the sub-bandgap emission, coupled with bandgap 
peaks, indicates thatradiative recombination combines 
with indirect intraband transitions (Bremsstrahlung) to 
produce the observed photon spectrum. The observed 
polarization of the light from GaAs/AlGaAs field- 
effect devices (MESFETs) may be attributed to a 
dominance of forward-directed electron scattering, an 
effect predicted theoretically. The spectrum of 
radiation fromthe 2DEG of GaAs/AlGaAs HEMTs 
explicitly shows the energy-dependent real space 
transfer effect, and the “vertical” wansfer of 2D 
conduction electrons to other device layers may be 
studied in detail. 
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These measurements, presently used to study 
transport in 2D systems, are currently being extended 
to examination of high-field transport in 1D devices. 
The electroluminescence spectra will likely provide a 
useful means to investigate the high-field and non- 
equilibrium transport in such lower-dimensional 
structures. 


1-4 High Injection Effects in GaAs/AlGaAs 
Quantum Wells: Spontaneous 
Recombination and Band-Gap 
Renormalization, N. Kirstaedter, E. H. Béttcher, 
M. Grundmann, and D. Bimberg, Institut far 
Festkérperphysik der Technischen Universitat, Berlin, 
West Germany. A comprehensive investigation of 
the charge carrier recombination and the band-gap 
renormalization (BGR) in GaAs/GaAlAs. single 
quantum well (SQW) structures is presented. The 
results are of particular importance for the modeling 
of QW-based semiconductor devices. By, measuring 
the differential lifetime, the spontaneous emission 
rate and the chemical potential as functions of the 
carrier injection rate in GaAs/AlGaAs SQW GRIN- 
SCH lasers, the carrier density dependence of the 
recombination rate and the BGR is determined in a 
sheet density range of 3x10!! cm-2 up to 7x1012 
cm~2, The radiative and non-radiative components 
contributing to the total emission rate are identified. 
It is shown that the radiative band-to-band 
recombination coefficient decreases significantly with 
increasing density. At low densities, the non- 
radiative contribution is dominated by interface 
recombination. At high densities, Auger- 
recombination becomes an important factor. 

For a typical 7 nm GaAs/Gag2Ali;gAs QW 
structure at room temperature, the BGR increases 
from -20 meV to -70 meV. In contrast to previous 
investigations, which are based upon the fitting of 
photoluminescence or gain spectra, the number of 
fitting parameters for the extraction of BGR data is 
significantly reduced. In addition, it is shown that an 
appropriate band model for the QW, which takes into 
account band mixing effects of the valence bands, has 
to be used for the calculation of the electron and hole 
quasi-Fermi potentials, 


1-5 Novel I-V Characteristic of 
Superlattices, Ping Ao, University of Illinois, 
Urbana, IL, USA and Jgrgen Rammer, Institut for 
Fysikk, Norges Tekniske Hogskole, Universitetet i 
Trondheim, Trondheim, Norway. In artificial 
Structures, such as superlattices, arbitrary 
relationships between bandstructure parameters can be 
achieved, enforcing the need for a reinvestigation of 
interband transitions of a crystal electron originally 
studied by Zener. We report the existence of non- 
exponential corrections to the Zener tunneling rate 
and demonstrate their importance in parameter 


regimes relevant fox quantum transport properties of 
superlattices such as Bloch oscillations. Novel 
oscillations in the I-V characteristic of a superlattice 
are predicted with periods determined by the band 
structure parameters. 


1-6 Monte Carlo Simulation of Lateral 
Surface Superlattices in a Magnetic Field, 
T. Yamada, A. M. Kriman, and D. K. Ferry, Arizona 
State University, Tempe, AZ, USA. When a 
low magnetic field is applied to an electron gas in a 
crystal, the spectrum consists of sharp, closely spaced 
Landau levels, With increasing field, the finiteness of 
the lattice constant leads to nonparabolicity 
corrections, which broaden and shift these levels. 
When the cyclotron radius 7, becomes comparable to 
the lattice constant a, a perturbative approach fails. 
Hofstadter! and others have studied the full range of 
behaviors that occurs for general values of r¢/a. Fora 
square lattice in the strong !attice-potential limit, the 
energy spectrum plotted as a function of the number 
of flux quanta per ¢ unit vell is fractal, with a period 
of unity. A single period of this spectrum takes the 
form of a “butterfly.” There are well-defined energy 
gaps, but the energy bands are extremely sensitive, 
nonanalytic functions of the field. 

In real crystals, the lattice constants are so small 
that the magnetic fields corresponding to the above 
condition are unachievably large - on the order of 
108G, However, this condition can now be realized 
in two-dimensional (2D) electron gases with 
artificially engineered lattices. In these lateral surface 
superlattices (LSSLs), 2D superlattice potentials can 
be achieved with periods on the order of 0.lum, 
leading to the 10° reduction in the field required. 
LSSLs constructed on MESFET and HEMT 
structures with meshed gate electrodes have confirmed 
the periodicity in magnetic field, and displayed a 
variety of other transport properties.» 

In order to study these structures, we perform a 
semiclassical Monte Carlo simulation of quasi-2D 
electrons described by a Stern-Howard wave 
function,* including electron-electron interaction, and 
scattering by acoustic phonons, optical phonons, and 
impurities. The electron-electron interaction is treated 
using a molecular dynamics technique.” We model 
GaAs-based LSSLs at 4.2K, with electron 
concentrations of 1.4x10!9 cm-?, and a superlattice 
potential described by Vo[cos(2nx/a)+cos(2ny/a)] 
with Vo= 10meV and superlattice period a = 0.16pm. 

The diffusivity is evaluated as a function of the 
magnetic field. It changes abruptly with small 
changes in the magnetic field, suggesting a fractal 
structure. The fractal structure found in the quantum 
case can be understood in terms of commensurability 
of two length scales - the superlattice potential period 
and the cyclotron radius. In the classical case studied 


Te oe 


here, a commensurability issue is still present which 
concerns the cyclotron frequency and the frequency of 
oscillation in the superlattice-potential minimum. 
1D, Hofstadter, Phys. Rev. B 14, 2239 (1978). 

2 G. Bernstein and D. K. Ferry, Z. Phys. B 67, 449 
(1987). 

3 R. A. Puechner, J. Ma,.R. Mezenner, W..-P. Liu, 
A, M. Kriman, G. N. Maracas and D. K. Ferry, Surf. 
Sci, 228, 520 (1990). 

4 F. Stern and W. E. Howard, Phys. Rev. 163, 816 
(1967). 

> P, Lugli and D. K. Ferry, Phys. Rev. Lett. 46, 
594 (1985). 


1-7 Far-Infrared Photoconductive Response 
of a One-Dimensional GaAs-(Ga,AlAs 
Channel in Zero and High Magnetic Fields, 
T. J. B. M. Janssen,2 N. K. Patel,! J. Singleton,? 
M. Pepper,! H. Ahmed,! D. G. Hasko,! R. J. 
Brown,! J, A. A. J. Perenboom,? G. A. C. Jones,! 
J. E. F. Frost,! D. C. Peacock! and D. A. Ritchiel, 
(University of Cambridge, Cambridge, United 
Kingdom; University of Nijmegen, The 
Netherlands). We have studied the photoconductive 
response of a one dimensional (1D) channel to 
monochromatic far-infrared (FIR) radiation from a 
laser (wavelengths 100-1200 jim) as a function of 
applied magnetic field at temperatures ~ 0.4 K. The 
channel is defined by a negative voltage applied to 
two gates 0.3 pm wide and 0.3 pum apart on top of a 
high mobility GaAs-(Ga,Al)As heterojunction: the 
transport through the channel is ballistic, and in zero 
magnetic field the resistance is quantised to a value 
hi2ie2, where i is the number of occupied 1D 
subbands. In a magnetic fielld applied perpendicular 
to the plane of the heterojunction, the resistance as a 
function of gate voltage exhibits plateaux due to 
conductionthrough the channel by a finite number of 
edge states, 

With a constant current through the channel, and in 
arbitrary magnetic field, the photovoltage across it 
due to the FIR exhibits oscillations as the gate 
voltage, and hence the number of 1D subbands or 
edge states within the channel, is varied: the 
oscillations are at positions related to the resistance 
plateaux, with the peakvalue of the photovoltage 
occurring to the low conductivity side of each plateau. 
It appears this non-resonant response results from an 
alteration of the electron energy distribution on cither 
side of the channel due to the heating effect of the 
FIR, 

At certain values of the magnetic field and gate 
voltage, however, the channel shows a strong 
resonant photoresponse, which is thought to 
consist of two components: 

i) when the FIR is resonant with the Landau level 
spacing outside the channel, the first empty Landau 


level above the Fermi energy my be populated by 
photoexcited electrons: these photoexcited electrons 
will have access to extra states for conduction through 
the channel; . 2 
ii) when the FIR is resonant with the edge state 
spacing within the channel, electrons may be excited 
to states above the Fermi energy in the channel. 
These results may be teated within the Biittiker 
formalism for phase coherent conductance between 
different electron reservoirs: within the reservoirs the 
electron energy distribution is altered by the 
absorption of FIR radiation. The implications of 
these observations for ieliftting 9 eevon in the 
Landau levels and channel states above the Fermi 
energy will be discussed. - 


1-8 Contribution of Ballistic Electrons in 
Vertically Integrated Resonant Tunneling 
Diodes, R. E. Carnahan, J. J. L. Rascol, K. P. 
Martin, and R. J. Higgins, Georgia Institute of 
Technology, Atlanta, Georgia, U.S.A., L. A. Cury 
and J. C. Portal, LPS-Institut National des Sciences 
Appliquees, Toulouse cedex and Service National des 
Champs Intenses-Centre National de la Recherche 
Scientifique, Grenoble cedex, France, B. G. Park, E. 
Wolak, K. L. Lear, and J. S. Harris, Jr., Stanford 
University, Stanford, California. We have made a 
systematic study of the current-voltage (I-V) 
characteristics of vertically integrated resonant 
tunneling diodes (RTDs) in the presence of a 
transverse magnetic field (B L J) up to 18 T and at 
temperatures below 4.2 K. The unit RTD consisted of 
17 A AlAs barriers and a 45 A GaAs well, surrounded 
on both sides by a 100 A undoped GaAs spacer layer. 
By varying the n-doped GaAs separation layer 
between two identical RTDs, we were ahle to 
systematically study the effect of ballistic electrons 
on the device. In the decoupled case (1000 A doped 
separation layer), relaxation of the electrons resulted 
in two negative differential resistance (NDR) peaks 
which both behave as single RTD peaks under the 
influence of B, However, the non-negligible 
minority ballistic electrons in the intermediate doped 
spacer layer (500 A) sample had their incident kinetic 
energy magnetically tuned to be resonant with the 
second quantized level in the collector RTD, causing a 
repartition of the potential across the RTD. Thus the 
second NDR transition (by the majority thermalized 
electrons in the current) was shifted to a higher bias 
by 400 mV. The third sample (with no doped 
separation layer) was dominated by ballistic transport 
and showed interference features for biases below the 
first peak voltage. The transverse magnetic field 
enabled the first peak current to dramatically increase 
with a slight increase in bias. An anomalous strong 
decrease in both the second peak voltage and current 
was ohserved. Additionally, a feature we associate 
with tunneling through the X valley in the AlAs 


collector. RTD barriers is observed for B > 10 T. We 
have shown that, even if the major portion of the 
current is thermalized, strong effects due to ballistic 
électrons are observed in integrated heterostructures. 


1-9 Elastic and Inelastic Resonant 
Tunneling in an- Imperfect Superlattice, P. 
Hyldgaard and A. P. Jauho, H. C. Orsted Institute, 
Copenhagen, Denmark. A model calculation of 
vertical transport through an imperfect superlattice is 
presented. A tight-binding model is used to describe 
transport in the superlattice (SL). The imperfection 
we study can be modeled by a double barrier. We 
include effects of (dispersionless) longitudinal optical 
phonon scattering in the quantum well. 

The model calculation includes the miniband that 
arises in an SL, and resonant tunneling within the 
miniband is found. Phonon-assisted resonant 
tunneling is found to:be operative, and a phonon 
satellite is identified. 

The calculations are done for experimentally 
realizable parameters, and constitute predictions for 
the experimental situation of two ideal superlattices 
interrupted by a set of neighboring higher/thicker 
barriers, with the entire imperfect superlattice situated 
in the base of a transistor. 


1-10 Magnetic Field Effect on Resonant 
Tunneling in Short-Period Superlattices, A. 
M., Vasilev and I. N. Uraltsev, A. F. Joffe Physico- 
Technical Institute, Leningrad, USSR. There has 
recently been a dramatic increase in research activities 
in the transport of carriers perpendicular to the layers 
of the strongly coupled superlattices (SLs), which has 
been shown to occur through extended, Bloch-type 
states.!:2. The perpendicular transport is expected to 
be suppressed in the presence of magnetic field 
parallel to the SL layers due to the wave function 
squeezing resulting in carrier localization when the 
magnetic length becomes less than the SL period. 
Destroying of the SL extended states has been 
predicted to occur at remarkably low magnetic ficlds 
as a result of the untuning of the resonance between 
electron states when cyclotron energy exceeds the 
miniband halfwidth.> 

We report the observation of the magnetic-field- 
induced suppression of the heavy-hole and electron 
transport through the extended states in 
GaAs/Al,Ga}_xAs SLs, A magnetic field applicd 
along the SL layers is found to produce a remarkable 
effect on the SL luminescence used for study of 
ambipolar transport of photocarriers in SLs with an 
enlarged QW. The efect is shown to depend 
dramatically on the SL period, d, for d<50A to result 
from the magnetic-field-induced transition from 
coherent Bloch transport to the transport through 
localized states of the heavy holes, To demonstrate 
destroying of the resonance between clectron states in 
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the presence of magnetic field we have measured the 
diamagnetic shift of the heavy-hole exciton in the 
SLs with narrow electronic minibands. In SLs with 
d<100A we have observed the strong quenching of the 
excitonic shift associated with the magnetic-field- 
induced localization of electrons. 

1B. Deveaud, J. Shah, T. C. Damen. B. Lamben, A. 
Chomette, and A. Regreny, JEEE J. Quant. Elec. 24, 
1641 (1988). 

2 P. S. Kop’ev, R. A. Suris, I. N. Uraltsev, and A. 
M. Vasiliev, Solid State Commun. 72, 401 (1989). 
3 A. M. Berezhkovskii, and R. A. Suris, Sov. Phys. 
JETP 86, 109 (1984). 
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2-1 Thermopower in Scanning Tunneling 
Microscope Experiments, J. A. Stévneng, 
Institutt for fysikk, Norges Tekniske H¢gskole, 
Universitetet i Trondheim, Trondheim, Norway, and 
P, Lipavsky, Ohio State University, Columbia, 
Ohio, U.S.A. We present a theory for the 
thermopower observed in scanning tunneling 
microscopy (STM), The lateral variation of the 
thermopower is found to depend on the logarithmic 
derivative of the local sample density of states at the 
Fermi level. We also derive a relation between the 
thermopower and the nonlinear conductance, The heat 
transfer dueto the tunneling electrons obeys the 
Weidemann-Franz law. 


2-2 Effectiveemass Boundary Conditions 
for Strained Heterostructures, G. T. Einevoll 
and P. C. Hemmer, Jnstitutt for Fysikk, Norges 
Tekniske Hogskole, Universitetet i Trondheim, 
Trondheim, Norway. The question of the appropriate 


form and validity of the effective-mass approximation: 


when applied to strained heterostructures is addressed. 
Complications arise from the position-dependence of 
the effective mass and lattice constant, and the form 
of the appropriate kinetic operator in the effective- 
mass Hamiltonian is not a priori given. For the one- 
band theory we propose the following two-parameter 
family of hermitian kinetic operators 


2 
Hkin = - Emad VmBba-25Vm%a5 , (1) 


where m and a are the, in general, position-dependent 
effective mass and lattice constant respectively, and 
one has the obvious requirement 20+ =-1. The 
corresponding effective-mass boundary conditions to 
use at abrupt heterointerfaces are uniquely determined 
by the form of Hgjn, and it follows that the quantities 


m%a5o and mb+%a-5Vo, where 6 is the effective 
mass wave function, must be conserved at. material 
interfaces, In our approach the parameters o, 8, and y 


re 


are determined by comparing effective-mass results 
with exact results for simple-microscopic models for 
abrupt: heterostructures. Within our calculational 
scheme we determine uniquely @ = 0 and B =-1 

while two distinct values, 0 and -1, are observed for 
5. Based on qualitative similarities between the test 
models and realistic systems, we propose for 
calculations on conduction-band states the boundary 
conditions listed above with a = 0, 8 =—1 and 6 =— 


1, When a one-band effective-mass theory is- 


applicable for valence-band states, however, the 
parameter values « = 0, 8 = —-1 and 5 = 0 should be 
used. “i N 

2-3 Berry’s Phase and Persistent Charge 
and Spin Currents in Textured Mesoscopic 
Rings, Daniel Loss, Paul Goldbart, and A. V. 
Balatsky, University of Illinois, Urbana; Illinois, 
U.S.A. The quantum orbital motion of eléctrons in 
mesoscopic normal metal rings threaded by a 
magnetic flux produces striking interference 
phenomena such as persistent currents and the 
Aharonov-Bohm effect. Sin::larly, when-a quantum 
spin adiabatically follows a magnetic field which 
rotates slowly in time, the phase of its state vector 


acquires an additional contibution known as the Bony : 


phase. 

The purpose of this paper is to explore the 
combination of these two quantum phenomena by 
examining the interplay between orbital and spin 
degrees of freedom for an electron moving in a 
mesoscopic ring. To this end, we consider a ring 
which is placed in a classical, static, inhomogeneous 
magnetic field, i.e., a texture. Zeeman coupling 
between the electron spin and this texture results in a 
Berry phase and, as a consequence, the system 
Supports persistent equilibrium spin and charge 
currents, even in the absence of conventional 
electromagnetic flux through the ring. We mention 
the possibility of analonous persistent mass and spin 
currents in normal 3He and spin-polarized H. 


2-4 A Numerical Method for Ihe 
Calculation of Transient Response in 
Mesoscopic Devices, Leonard F. Register and 
Umberto Ravaioli, University of Illinois, Urbana, 
Illinois, U.S.A, We present a numerical method for 
the solution of the 2-D time-dependent Schrocdinger 
equation, suitable for the investigation of transients 
in mesoscopic devices under ballistic conditions. The 
features of the method are: (1) a tight-binding 
formulation of the quantum mechanical Hamiltonian; 
(2) 2-D absorbing boundary conditions to simulate 
the discretized open system; (3) Crank-Nicholson 
scheme for the evaluation of the space dependent 
Operator, solved with an Alternate Direction Implicit 
(ADI) strategy. We believe that this is the first 
application of 2-D absorbing boundary conditions for 


the Schroedingér equation. The ADI scheme decouples 
the discretized linear system of equations into subsets 
requiring only the solution of tridiagonal systems. In 
addition, since all dependencies of data are eliminated, 
the code can be successfully vectorized [1]. The 
method allows us to investigate separately the various 
energy components, instead of having to resort to 
wavepackets, and this results in better flexibility for 
the future application and comparison of different 
formulations for the injecting conditions at the 
contacts. To illustrate the features of the algorithm, 
we present results of switching calculations for an 
ideal T-structure proposed earlier as a quantum 
modulated transistor. Although the structure is 
modelled as an ideal electron waveguide, the essential 
features of the quantum interference effects during the 
tum-on transient can be appreciated and understood. 

(1) J. R. Barker, J. Pepin, M. Finch and M. 
Laughton, Solid State Elec., 32, 1155 (1989). 

(2] F. Sols, M. Macucci, U. Ravaioli and K. Hess, J. 
Appl. Phys., 66, 3892 (1989): 


2-5 Dephasing by a Dynamic Environment, 
D. Loss and K. Mullen, University of Illinois, 
Urbana, Illinois, U.S.A. We investigate the manner 
in which quantum interference is suppressed whena 
particle interacts with a localized, dynamical 
environment. To do so we examine a model with 
two classical paths along which an electron can 
travel, and allow it to interact with a bath of 
harmonic oscillators on one path, and travel freely on 
the other, In particular we show that the quantum 
fluctuations of the path of the particle can couple to 
the environment and thus lead to dephasing, and 
calculate the dephasing time in the high temperature 
limit. We compare this result to other views of how 
propagating electrons lose phase coherence. 


2-6 Scattering in Nearly-Clean Mesoscopic 
Structures, A. M. Kriman, B. S. Haukness and D. 
K. Ferry, Arizona State University, U.S.A. We 
study theoretically the effect of small numbers of 
elastic scatterers on ballistic transport in low- 
dimensional microstructures. Both numerical and 
analytical methods are used, within a transfer matrix 
formulation. With lateral confinement, delta-function 
and other extremely sharp models of a single defect 
lead to sharp resonances when such defects are well 
isolated, occurring when total electron energy equals 
any miniband energy, Unit transmission probability 
is approached when the scattering defect is small, and 
far from other scatterers, even if the scatterer is strong 
enough to decrease significantly the conductance away 
from resonance. Resonances occur for all shapes of 
confining potential. 

At low energies, single-barrier structures with a 
single nearby defect exhibit resonant transmission 
similar to that of double-barrier resonant tunneling 


diodes (DBRTDs), with an approximate scaling 
behavior that relates transmission for defects at 
different distances to that at a fixed distance with 
different energy scales. In ordinary DBRTDs, the 
position of the transmission peak is strongly affected 
only by defects lying within the quantum well region. 
The height of the transmission peak is very sensitive 
to the positions of defects within that region, which 
act essentially as a probe of the resonance wave 
function. Defects in front of a DBRTD also affect the 
valley current by modifying the longitudinal 
component of the incident momentum. 


2-7 The Effect of Dissipation on Phase 
Periodicity and the Quantum Dynamics of 
Josephson Junctions, D. Loss and K. Mullen, 
University of Illinois, Urbana, Illinois, U.S.A. We 
consider systems described by compact, periodic 
variables, such as Joseplison junctions and Small 
normal metal rings. We examine how they can be 
coupled to a heat bath, via either momentum or 
position coupling, and how the two descriptions can 
be related by a unitary transformation. We show how 
it is critical to wansform not only the Hamiltonian, 
but also the initial conditions. We then demonstrate 
that for certain types of initial conditions, paths of 
different winding number can interfere. Still other, 
ostensibly reasonable initial conditions, have no such 
interference. We conclude with a discussion of 
appropriate models to describe periodic systems. 


2-8 The I-V Characteristic of a Resistively 
Shunted, Small Capacitance Josephson 
Junction, V. Bubanja, A. Maassen van den Brink, 
D. V. Averin® and G. Schén, Delft University of 
Technology, Delft, The Netherlands. We consider the 
low temperature behaviour of a small capacitance 
superconducting junction which is shunted by an 
ohmic shunt. The dynamics of such junctions was 
Studied in ref. [1] within the framework of a 
deterministic single-band model. We allow finitic 
values of Zener tunneling probability and include 
transitions bewtween different quantum states caused 
by the shunt. A similar situation, but where the 
dissipation is due to quasiparticle tunneling, has been 
discussed in refs. (2] and [3]. Here we concentrate on 
the limit of negligible quasiparticle tunneling, For 
currents less than e/RsC (Rg denotes shunt resistance, 
C junction capacitance) all the current is going 
through the shunt, so that the I-V curve is linear. 
Beyond it follows a branch with negative differential 
resistance due to Bragg reflection and due to the 
Cooper pair tunneling induced by the shunt. For 
large currents Zener tunneling prevails and the I-V 
curve bends over to the ohmic asymptote. 

permanent address: Moscow State University, 
Moscow, U.S.S.R. 


(1] K. K. Likharev and A. B. Zorin, J. Low Temp. 
Phys. 59, 347 (1985). 

{2] D. ¥. Averin and K. K. Likharev, J. Low Temp. 
Phys. 62, 345 (1986) and to be published 

(3) U. Geigenmiiller and G. Schén, Physica B 152, 
186 (1988). 


2-9 Single-Electron Tunneling in Point- 
Contact Tunnel Junctions, R. T. M. Smokers, 
P. J. M. van Bentum and H. van.Kempen, University 
of Nijmegen, Toernooiveld, Nijmegen, The 
Netherlands, Charging effects like Single-Electron 
Tunneling or the Coulomb staircase can be observed 
in point-contact tunnel junctions on a wide range of 
materials, Thsi indicates that the charging energy 
Q2/2e can be the dominant term in the Hamiltonian 
describing the tunneling in these junctions. Recent 
theoretical developments have elucidated that this can 
only be the case if the high-frequency impedance of 
the junction environment is large enough to 
effectively decouple the junction from its long-range 
electromagnetic environment (e.g. stray capacitances). 
For a series connection of tunnel junctions, as in the 
case of tunneling through small isolated paticles, this 
condition is obviously satisfied. The observations in 
single junctions however cannot be readily 
understood. To investigate the influence of the 
electromagnetic environment we have performed 
tunneling measurements on surface-doped Si, at 
temperatures close to the metal-insulator transition. 
We observe clear S.E.T.-behavior, with capacitances 
of order 10~!7F, but find no appreciable effect when 
changing the series resistance of a point contact from 
=50kQ to=80MQ. Decoupling from stray 
capacitances may be due to a relatively slow hopping 
conduction through impurity states in the direct 
vicinity of the point contact. 


2-10 A Generalized Impact-Ionization 
Model for High-Energy Electron Transport 
in Si with Monte Carlo Simulation, 
Rossella Brunetti, Universita di Modena, Italy. A 
new model for impact ionization in Si is presented, 
which goes beyond the limitations of the Keldysh 
forrnula and is based on a more realistic scheme 
developed starting from a first-order perturbation 
theory. This scattering mechanism is modeled by an 
extended band structure which includes many bands 
for clectrons and one band for holes in a finite 
Brillouin zone. Some processes have been identified 
to bring the dominant contribution to the scattering 
probability, in the present approach, foir clectron 
energies up to 3eV. Expressions for the differential 
and integrated scattering probabilities have been 
obtained which are consistent with the band model 
and can be included in a Monte Carlo simulation of 
the electron gas. Results for transport quantities are 
presented for a bulk material in presence of 
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chomogeneous and static electric fields under physical 
conditions where impact ionization influences the 
-carrier dynamics. A comparison with theoretical and 
experimental data from the literature is also given. 


2-11 Ambipolar Perpendicular Transport in 
a Semiconductor Slab: Influence of Optical 
Phonons, K. Scheller, T. Held, and G. Mahler, 
Universitat Stuttgart, West Germany. We investigate 
the perpendicular transport of an optically excited 
carrie plasma in.a-thin semiconductor slab. The 
excitation is done by a monochromatic laser, 
homogeneously illuminating one surface. A kinetic 
description by means of the Boltzmann equation is 
used. On length scales greater than the absorption 
length and on time scaies greater than the carrier- 
carrier scattering time the distribution function is 
approximated by heated displaced Maxwellians. We 
focus on the role of optical phonons for the transport 
and relaxation properties of the plasma. 


2-12 Phonon scattering and energy 
relaxation in two, one and = zero 
dimensional electron gases, U. Bockelmann 
and G. Bastard, Laboratoire de Physique de la Matiére 
Condensée de l'ENS Paris, France. We report on 
calculations of intra- and inter-subband phonon 
scattering in quantum confined electron gases based 
on lattice matched InGaAs/InP quantum wells. The 
emission of longitudinal acoustic (LA) phonons from 
two, one and zero dimensional electron gases are 
compared by studying the scattering times as a 
function of the width of the lateral confinement. 
Longitudinal optical (LO) phonon scattering in 
quantum wells and wires are discussed using a 
phenomenological broadening of the 1 D density of 
states. The energy relaxation rates of heated electron 
gases due to phonon emission and absorption are 
calculated for lattice temperatures Tj between 0.3 and 
20K as a function of the electron density. In the 1 D 
systems the characteristic quantities (scattering rates 
+71, mean LA phonon energies (Eph), energy 
relaxation rates P) exhibit oscillations around their 
corresponding 2D values, which reflect the 1 D 
density of states. On the other hand, the LA phonon 
scattering in OD systems becomes increasingly 
quenched with increasing quantization energies. 
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homogencous and static electric fields under physical 
conditions where impact ionization influences the 
carrier dynamics. A comparison with theoretical and 
expcrimental data from the‘literatuié is also given. 


2-11 Ambipolar Perpendicular Transport in 
a Semiconductor Slab: Influence of Optical 
Phonons, K. Scheller, T. Held,-and G. Mahter, 
Universitat Stuttgart, West Germany. We investigate 
the perpendicular transport of aa optically excited 
carrie plasma in a thin semiconductor slab. The 
excitation is done by a monochromatic laser, 
homogeneously illuminating one surface. A kinetic 
description by means of the Boltzmann equation is 
used. On length scales greater than the absorption 
length and on time scales greater than the carrier- 
carrier scattering time the distribution function is 
approximated by heated displaced Maxwellians. We 
focus on the role of optical phonons for the transport 
and relaxation properties of the plasma. 


2-12. Phonon scattering and _ energy 
relaxation. in two, one and = zero 
dimensional electron gases, U. Bockelmarn 
and G. Bastard, Laboratoire de Physique de la Matiére 
Condensée de l'ENS Paris, France. We report on 
calculations of intra- and inter-subband phonon 
scattering in quantum confined electron gases based 
on lattice matched InGaAs/InP quantum wells. The 
emission of longitudinal acoustic (LA) phonons from 
two, one and zero dimensional electron gases are 
compared by studying the: scattering times as a 
function of the width of the lateral confinement. 
Longitudinal optical (LO) phonon scattering in 
quantum wells and wires are discussed using a 
phenomenological broadening of the 1 D density of 
slates. The energy relaxation rates of heated electron 
gases duc to phonon emission and absorption are 
calculated for lattice temperatures Tj between 0.3 and 
20K as a function of the electron density. In the 1 D 
systems the characteristic quantities (scattering rates 
tl, mean LA phonon energies (Eph), energy 
relaxation rates P) exhibit oscillations around their 
corresponding 2D values, which reflect the 1 D 
density of states. On the other hand, the LA phonon 
scalicring in OD systems becomes increasingly 
quenched with increasing quantization energies. 


2-13 Reduction of the conductance 
fluctuation amplitude in mesoscopic 
GaAs/GaAlAs heterojunctions due to spin- 
orbit scattering, M. W. Keller, O. Millo, S. J. 
Klepper, S. Xiong, A. D. Stone, and D. E. Prober, 
Department of Applied Physics, Yale University, 
New Haven, CN 06520, and R. N. Sacks, United 
Technologies Research Center, East Harford, CN 
06108. We have studied week-localization and 
conductance fluctuations in mesoscopic two- 


dimensional GaAs/GaAIAs heterojunctions for T=7 K 
to60 mK. The weak localization data show that.the 
spin-orbit scattering rate exceeds the inelastic rate 
below ~ 2K. At the same time, we find a significant 
reductionzin the conductance fluctuation amplitude- 
below ~ 2K, as compared to that extrapolated from 


~ the high temperature region, where the expected T-1/2 


behavior is observed. Our results agree well with 
calculations for-the effect of spin-orbit scattering on 
the universal conductance fluctuation amplitude. The 
effect of spin-orbit scattering on the magnetic 
correlation range is also discussed. 


2-14 Sensitivity of conductance 
fluctuations in mesoscopic devices to 
individual elastic scatterers, S.J. Klepper, O. 
Millot, M. W. Keller, and D. E. Prober, Department 
of Applied Physics, Yale University, New Haven, 
CN 06520": and R. N. Sacks, United Technologies 
Research Center, East Harford, CN 06108. We have 
added elastic scatterers in a controlled fashion to 
mesoscopic wires fabricated-from an AlGaAs/GaAs 
heterojunction, through the photoionization of DX 
centers in the AlGaAs. These devices are in the 
diffusive regime. The attached leads are optically 
masked, so that the Fermi parameters of the 2D 
electron gas do not change during infrared (IR) 
illumination. We are able to resolve switching 
events in the device conductance due to the addition of 
single scatterers. By studying the evolution of the 
conductance fluctuation characterisitic "magneto- 
fingerprint" of a sample between successive IR 
illuminations, we are also able to observe the 
Statistically-averaged effect on device conductance of 
adding a number of scattering sites. From this we 
can determine the rms contribution of a singlc elastic 
scatterer. Our results for the amplitude and 
termperature dependence of these conductance changes 
are consistent with theoretical predictions.! 

* Supported by NSF DMR 8505539 

* Weizmann Fellow 

1. S. Feng et al., Phys. Rev. Lett. 56, 1960 (1986) 
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_ AN INTRODUCTION 
TO CHARGE QUANTUM TRANSPORT IN SEMICONDUCTORS 
AND NUMERICAL APPROACHES 


Fausto Rossi, Rossella Brunetti, and Carlo Jacoboni 
Dipartimento di Fisica dell’ Universita, Via Campi 213/A, I-41100 Modena, Italy 


INTRODUCTION TO QUANTUM TRANSPORT 


Since several years it has been recognized that the features of the new phys- 
ical systems provided by present-day technology require a quantum treatment of 
electron transport. For recent reviews on the subject see Grubin et. al., 1988 and 
Ferry et al., 1990. 

We shall not repeat here the reasons for such a belief. We shall rather con- 
centrate on an introduction to the basic ideas of quantum transport theory and 
to a brief account of some of the efforts that are presently being made to obtain 
numerical results for an estimate of the main quantum effects or, at least, for 
a test of new numerical approaches. In fact, even though the main theoretical 
concepts and methods of quantum transport have been developed long ago, they 
hardly yielded any numerical result to compare with experimental data, owing to 
the complexity of the mathematics involved and to the lack of clear experimental 
evidences. 

In the theory of quantum transport, many concepts and techniques are used, 
such as density matrix, Green functions, path integrals and Wigner function, which 
may be not familiar to the non-experts, and for this reason we shall try to give in 
this section a quick sunimary of such concepts. In particular, we shall enphasize 
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their differences, similarities, and mutual relationships. In the following sections | 
these concepts will be somewhat developed and some numerical applications ea 


be presented, without any attempt of completeness. 


The quantum description of a physical system is given, in the Schrodinger 
picture, by the state vector 


|®(z)) ~~“ (1) 

as a function of time ¢. 
If the problem is not completely specified, we have to use the concepts of 
statistical physics to deal with our incomplete knowledge, and the mathematical 
instrument in this case is the density-matrix operator (see, for example Ter Haar, 
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1961): 
p(t) = |B(t)){9(¢)| (2) 


where the overbar indicates an average‘to be performed over a suitable statistical 
ensemble that accounts for our partial knowledge of the system. 


In a given basis |¢,) the density-matrix operator has the following matrix 
elements 


Pum = ($n|®)(]dm) = cnet, (3) 

where c, are the coefficients of |®) in the given basis. 
The state vectoz or the density matrix can be used, for a pure state or an 
_ ensemble, respectively, to obtain the average result of a measurement of a quantity 


A: 


(A) = (8].4]S) (4) 
(A) = (@[AJS) = Sock Anmem = Tr(pA) . (5) 


If the system under examination can be devided into a “small” part of interest 
with coordinates z and a “large” part formed by an external interacting system 
with coordinates X, and if an observable A‘*) acts only on the variables z, the 
mean value of A‘*) is given by 


(A) > Tr(pA) = Yo (2, X, z', X)A)(2', 2) 


Poe, @ 
or 
where 
p(x,2') =) p(t, X,0',X) = Trx(p) (7) 


x 


is the reduced density matrix. 
The Wigner function (Wigner, 1932) is a special transform of the density 
matrix in the coordinate representation. In a one-dimensional case, for example, 


aN 


1 i 1 1 
fw(z,k,t) = on [eve *Y a(x = 9h” - 59) " (8) 


Therefore the Wigner function carries the same information as the density matrix. 


In order to solve a quantum problem with given initial conditions for the state 
vector or the density matrix, it is necessary to find the evolution operator U(¢,t;) 
from the initial time ¢; to the actual time ¢. In the Schrodinger picture, it evolves 
the state vector according to 


|®(¢)) = U(, &s)|B(E)) (9) 
2 
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and, therefore, the density matrix according to 


p(t) = U(t,ti)o(t) UN, t). (10) 


If {q} represents a complete set of compatible dynamical variables for-our system, 
Eq.(9) can be written in the {gq} representation as 


$(q,t) = / U(q,q',t, ti) B(g', ts)dq', (11) 


where-U(q,q',t,t;) is the matrix element of the evolution operator U(t,t;) in the 
basis {g}. 

- The fundamental dynamical equation is then a differential equation for the 
evolution operator: 


RLU (ist) = HU(t, ts), (12) 
with the initial condition 
U(ti, tz) = 1, (13) 


where H is the Hamiltonian of the system. 


The Feynman path-integral theory (see Feynman and Hibbs, 1965) is an al- 
ternative approach to quantum mechanics. It starts from the idea that all possible 
paths of a system from an initial state to a final one are to be considered as simul- 
taneously realized, and their amplitudes add up, rather than their probabilities as 
it would be in classical concepts, to give the probability amplitude of finding the 
system in the final state. The fundamental equation in this approach is an explicit 
expression for the evolution operator as an integral over all possible paths of the 
exponential of the classical action. If q represents a set of classical Lagrangian 
variables of the system, and g(r) one given trajectory, or path, from the initial 


values q; = q(t;) to the final valves gf = q(t), the evolution operator in Eq.(11) 
can be written as 


qtit . 
Ularsgiststs) =f Daf) eb Sat, (14) 
qirti 

Here Dq({r) indicates the integral over all paths that connect the initial state 


{qi,t;} to the final state {g7,2} and S[q(r)} is the classical action evaluated over 
each given trajectory in the integral: = 


State) = f° Har) a(r)s7)en (15) 


where L(q, 4,7) is the classical Lagrangian of the system. 

If a path-integral approach is followed and the integral in Eq.(14) evaluated, 
the resulting evolution operator can then be used to evaluate the evolution of a 
state wavefunction (as in Eq.(9) or (11)), or of the density matrix (as in Eq.(10)). 


As it regards the Green functions, they carry information on both the state 
of the system and its dynamical evolution (see, for example, Kadanoff and Baym, 
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1962; Rickayzen, 1980; Mahan, 1981). Several quantities-are often indicated with 
this name, and they are somewhat different, even though strictly related to each 
other. Some of them are given by an. avérage value of a commutator of field 
operators times a step function 9. More precisely, the single-particle retarded and 
advanced Green functions G(") and G‘*) are defined as 


Oe Aa 2!) = GnllWalre) Helo! t slate) (18) 
G(r, ¢,r',t') = = (Grrl, t), whe! t elbn)0(e =x t) (17) 


where |) is the state under consideration-in the Heisenberg picture, %,(r, t) is 
the field operator in the same picture, and @ is the Heaviside step function. The 
upper sign is to be used for Bosons and the lower sign for Fermions. An ensemble 
average must be included in case of statistical problems. 

In order to grasp the physical meaning of such functions let us recall that 
the expectation value of the product wl (r, t)Wo(r,t) at equal positions and times 
gives the intensity of the field in the state under consideration. The same product 
at different arguments gives, in the same way, the correlation. between the two 
amplitudes at different positions and times. If we look for the dynamical corre- 
lation, that is the propagator, without the information on the field intensity, we 
must subtract the product in reverse order. We may compare this result with the 
more familiar relations for harmonic-oscillator creation, annihilation, and number 
operators: 


ala=N, aalt=N-1, sothat ([a,at]=1 (18) 


Thus, the two Green functions above carry the information on which probability 
amplitude for the presence of a particle in r at ¢ corresponds to a unit amplitude 
in r' at ¢’. They are, therefore, the equivalent of a propagator, or of an evolution 
operator. They include the fact, however, that they are single-particle operators 
defined in a many-particle system, and therefore they contain a reduction of many 
degrees of freedom and include the effect, on each particle, of the interaction with 
all the other ones in the system. 

The average values of the single products of field operators discussed just now 
are themselves defined as other Green functions. More precisely, G7 and G< are 


GC? (oytyr' t!) = (Baler UC £9) (19) 
GS (rytyr' tt) = (Gh (re) ¥e( 1) (20) 
Clearly, the following relations hold: 
Gi) = a(t —2'\(G? -— GS)._. - (21) 
G) = ~9(' — 1)(G? — GS) (22) 


From what we have seen above it is understandable that G< is a correla- 
tion function of the field amplitudes at (r,t) and (r’,t'); thus, for equal times it 
corresponds to a reduced density matrix for a single particle. 
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‘A Wigner function can therefore-be defined starting from- fox (Mahan, 1981) 
which generalizes the definition given in Eq. (@) to a single-particle Wigner function 
for a many-body system: : 

Other Green functions can be defined for two or more particles. For such 
Green functions, including the single-particle ones, there are.no simple equations 
of motion since the propagation of a single particle depends on the entire many- 
body system. In fact it is possible to define a set of hierarchical equations of 
motion, where the equation for the single-particle Green function contains also 
the two-particle Green function; the equation for the two-particle Green function 
contains the three-particle Green function, and so on. 


o 
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THE MONTE CARLO DENSITY-MATRIX APPROACH 


The time evolution of the density matrix in Eq.(10) and the fundamental 
equation of motion (12) yield immediately the Liouville-von Neumann differential i 
equation for the evolution of the density matrix 


in = (2, pl. (23) 


As an example of application of the density-matrix approach we briefly sum- 
marize the work performed by the Modena group (Brunetti et al., 1989; Menziani 
et al., 1989; Rossi and Jacoboni, 1989) for the solution of the Liouville-von Neu- 
mann equation for the electronic density matrix in semiconductors. In principle, 
the method allows to evaluate the electronic density matrix as a function of time 
without any assumptions on the intensity and the duration of the electron-phonon 
interaction, as well as on the strength of the applied field. 

As starting point, let us consider a noninteracting electron gas in a semicon- 
ductor crystal, coupled to the phonon gas and to a constant and uniform electric 
field E. The system is assumed to be homogeneous, and its Hamiltonian is given 
by 

H=H.+He+ Hp + Hep (24) 


where H; is the term corresponding to an electron in a perfect crystal, Hz = eE-r 
is the term due to the electric field, and H, is the Hamiltonian of the free phonons; 
the electron-phonon interaction Hamiltonian H,, has the general forrir 


Hey = > ihF(q){age'"* — ate *4"} = Has + Hem (25) 
q 


where H,, and Hem refer to phonon absorption and emission respectively, and 
F(q) is a function of the phonon momentum q whose explicit form depends on 
the particular interaction mechanism. 

Let us consider the set of time-dependent basis vectors | Ko, {rq}, t) repre- 
sented by 


1; -i w{k(r 
penne he eM | tht (26 
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wheré k(t) =k. — si, and w(k) represents the electronic band structure. They are 
direct products-of electronic accelerated plane waves, or Houston waves-(Houston, 
1940), normalized to 1 over the crystal volume V,, and the phonon states | {ng}, t). 

The Liouville-von. Neumann equation that describes the time evolution of the 
density matrix p of the system in-the representation of the set in Eq.(26) contains 
only the perturbation Hamiltonian: : 


Pa) . . 
in o(2y2's#) = [Heps (tla, (27) 


where the symbolic compact notation z = (k,,{mq}) has been used. Eq.(27) has 
the same form we would have obtained in an interaction picture, even though 
the basis functions (26) are not eigenstates of the unperturbed hamiltonian. This 
is so because the Houston waves are solution of the time-dependent Schrodinger 
equation with the unperturbed Hamiltonian that -includes the electric field. 

If one is interested in the evaluation of expectation values of electron quantities 
which are diagonal in the electronic part-of the states in Eq.(26) only the diagonal 
elements p(z,¢) = p(z,z,t) of p must be determined. Furthermore, a diagonal 
initial condition for p decoupled in electron and phonon coordinates has been 
assumed. 

After a formal integration of Iiq.(27), a perturbative expansion for p is easily 
obtained by iterative substitutions: 


p(x,t) = p(z,0) +f dty (Hep, o()lti32,2 


t ty 
+f dt, / dt2[Hep, (Hep, p(O)|ea]t::2,2 Tee 
0 0 


= p(a2,t) + Ap Mx, 2) + ApM(z,t) +... 5 (28) 
where Hep = Hep: 


The trace over phonon coordinates can be performed exactly under the as- 
sumption that each phonon mode is involved only once during the transport pro- 
cess, which is equivalent to neglect hot-phonon effects. 

A diagrammatic representation of Eq.(28) allows us to regard each term of 
the perturbative expansion as a sequence of quantum processes which correspond 
to single scattering events in semiclassical transport. This constitutes “the starting 
point of the numerical quantum Monte Carlo (QMC) algorithm devised for the 
solution of the Liouville-von Neumann equation, which is based on random gener- 
ations of all possible sequences of processes associated to the different perturbative 
corrections, in the same way as the usual classical Monte Carlo (CMC) generates 
semi-classical scattering events. ~s 

Fig.l shows an example of a seventh order diagram, which includes also the 
electron-impurity interaction discussed below. It is worthwhile stressing that this 
QMC technique is based on a perturbative expansion of the density matrix ex- 
actly as the standard CMC technique is based on a perturbative expansion of the 
distribution function: a path with n scattering events in CMC is a term of the 
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n-th order in a_perturbative:expansion.of the Boltzmann equation in powers of the 
scattering rates (Poli et al., 1989). 

In what follows we present some typical results obtained with this QMC pro- 
cedure for different materials and physical conditions. The analysis has been 
devoted to very short times after the initial conditions, when quantum features 
are expected to be more relevant. This choice allows to include only few terms of 
the perturbation expansion given in Eq.(28), which, in turn, limits. the computer 
time to affordable values. 

Spécial emphasis has been given to the investigation of typical quantum ef- 
fects, such as “intracollisional field effect” (ICFE) and the “collisional broadening”, 
by comparison with the semiclassical case. 


A. Analysis of the electronic density matrix in presence of electric fields 


Numerical results have been obtained for the diagonal part of the electronic 
density matrix in presence of an arbitrary high electric field, starting from equi- 
librium conditions for the electron and phonon systems. Both electric field and 
electron-phonon coupling due to nonpolar optical phonons, are turned on at t=0. 
A simple-model semiconductor modeled on silicon with. a single spherical and 
parabolic band, and carriers interacting with one optical-phonon mode has been 
used. The working-conditions are T=20 K, B=150 kV/cm, t=50fs. 

Fig. 2 shows the quantum distribution function as a function of k along the 
field direction compared with the corresponding classical result obtained from the 
perturbative expansion of the Boltzmann equation by means of a “backward Monte 
Carlo” (BMC) procedure (Jacoboni et al., 1988). The distribution functions are 
peaked around the ballistic value. The particles that are scattered out of the 
ballistic trajectories are spread in a large volume of k space and cannot be seen 
in the figure. The quantum distribution is lower than the classical one. It can 
be shown that this correction is the net result of two larger quantum effects of 
opposite signs. In fact, if the effect of the field during the collision is neglected, 
we may obtain a big effect of non energy-conserving transitions, due to the short 
time considered. When the ICFE is added, it reduces the scattering efficiency by 
reducing the time of positive interference which occurs when the energy difference 
between initial and final states is equal to the phonon energy. 


B. Transient analysis of electron-impurity interaction 


The theoretical approach discussed before has been also extended to the anal- 
ysis of electron-impurity interaction, adding to the Hamiltonian in Eq.(24) an 
interaction term given by 


N 
He =>) o(e-1); (29) 


t=1 


where r;,r2...ry describe a given impurity configuration within the crystal vol- 
ume V, and y(r —r;) is the interaction potential due to the i-th impurity. 
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The matrix-elements between ko and kj of the electron-impurity Hamiltonian 
depend on the impuiity coordinates as 


N 


SS eWi(ko—ka)-i (30) 


i=1 


and in order to proceed in the calculation an ensemble average over all the impurity 
configurations {ri,r2...rv} is to be taken. These matrix elements have been 
incorporated into the QMC procedure described above. The diagonal part of 
the electronic. density matrix has been evaluated for the same simplified silicon 

_ model used previously, and the electron-impurity interaction is taken as a screened 
Coulomb potential. In Fig. 3 the quantum distribution function up to the fourth 
order is shown as a function of k along the field direction, compared with the 
corresponding classical result obtained from a BMC simulation (Jacoboni et al., 
1988). The distribution functions are peaked around the ballistic value. Here the 
quantum curve is higher than the classical one and this fact is again the net result 
of two larger quantum effects of opposite signs: on one hand, ICFE for a given 
transition (from k to k’) in momentum space tends to decrease the scattering 
efficiency; on the other hand, energy-non conserving transitions, allowed at so 
short times, may increase the number of the available final states. In this case, the 
prevailing effect is the first one since the particular form of the electron-impurity 
interaction itself favours low momentum-transfer and, therefore, energy-conserving 
transitions. 

Higher-order terms contain processes involving more than one scattering cen- 
ter (corresponding to several serniclassical scattering events), and processes involv- 
ing more than two vertices related to the same scattering center (corresponding 
to higher-order corrections to the Born approximation).. 


C. Quantum energy relaxation of photoexcited carriers 


The present QMC method has been also applied to the case of photoexcited 
electrons in bulk GaAs. The semiconductor model has been simplified to a sin- 
gle spherical and parabolic band similar to the one used before for Si and the 
interaction Hamiltonian includes only polar coupling to optical phonons. 

Electrons are generated at ¢ = 0 according to a distribution proportional to 
exp{—ale — €.|}, where ¢ is the electron energy, and q is an appropriate constant. 

For the sake of clarity we may consider first results of the classical theory. 
Fig. 4 shows the results obtained with CMC at ¢ = 100fs after excitation. The 
highest peak at 1000 K represents what is left of the initial distribution at ¢ = 0. 
Two secondary peaks are clearly seen, corresponding to electrons having emitted 
one or two optical phonons. Fig. 5 shows the corresponding result obtained with 
the present QMC (note the scale change). The initial distribution is diminished 
of a quantity very similar to that of the classical case. However electrons can be 
found, at t = 100fs, in a very wide range of energies, since energy needs not be 
conserved. The secondary peaks are not yet well formed. If we go towards longer 
times the secondary peaks appear also in the quantum result. 
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D. Quantunt analysis of drift velocity overshoot in GaAs 


Finally, the QMC procedure has been applied to study the drift-velocity over- 
shoot in GaAs. 

A simplified semiconductor model has been used, consisting of a I’ valley 
and four equivalent L valleys. The electron-phonon interaction is given by two 
scattering mechanisms: intra-valley polar optical and intervalley nonpolar optical 
phonons. The value. of the electric field has been taken as 40 kV/cm, and this 
choice is based on the analysis of the overshoot in classical terms. 

In Fig. 6 the drift velocity, as-given by a traditional CMC simulation, is been 
compared with its corresponding perturbative solution up to the second order in 
the scattering rates (i.e. the ensemble of trajectories with up to two scattering 
events). We can see that this perturbative solution gives us a good approximation 
of the drift velocity peak. Therefore we expect to be able to describe velocity 
overshoot in quantum terms with the perturbative expansion described above up 
to the fourth order in the hamiltonian. 

In fact, in the same figure the corresponding quantum drift velocity up to the 
forth order has been shown: it can be seen that under these conditions, quantum 
features are not very relevant. Because of energy-non-conserving transitions one 
could expect that the drift-velocity overshoot would be reduced in quantum theory, 
owing to an anticipation of the intervalley transitions. This does not occur however 
since the coupling constant is too weak. 


THE WIGNER-FUNCTION APPROACH 


In a sense the density matrix can be considered the quantum analogue of the 


classical distribution function, since its diagonal values in a fixed basis give the 
probabilities of finding a system of the ensemble in the corresponding eigenstate. 
However, the distribution function f(q,p) for a classical statistical ensemble gives 
us the probability density to find the system in the position (q,p) of the phase- 
space. Owing to the incompatibility of q and p, it is not possible to find a function 
with the same property in quantum mechanics. 

In this connection, we may rise the more general problem of defining a numeri- 
cal function A(q,p), where g and p are c-numbers, which is the classical analogue of 
a general quantum operator A(q,p), function of the operators q and p. A requisite 
could be that the sum of all possible values is the same for both quantities: 


Tr{a}= a | f Alap)dade, (31) 


This, however, is not enough to determine A since only an integrated property 
of it is given. On the other hand, if we require that the property expressed by 
Kq.(31) holds true also giving different weights to the different points in the phase- 
space, then the equation can be inverted. Here we have an arbitrariety in how to 
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choose these different weights. If Foutier-functions are chosén, we obtain.the Weyl 
correspondence rule. 
Let us then require that, for any real € and 77, 


Tr{ Aeiléa/ Atnp/A}y as ot i; / A(a, B)aise/* e08/h dad g, (32) 
Now the Fourier transform in the right hand side can be inverted: 


A(c,B).= x i e~ka/h ge / e~i08/* dy i dz i dz (x|Alz'\(z' [e%$9/*+7P/4] 12), 


, (33) 
By using 
eiléatnp) — einp/? gita ginp/2 (34) 
and 
ei?/*\n) = |x — 1), (35) 
we obtain, after streightforward calculations, 
A(a,p) = f e°/a—n/21Aq,e)la+n/2)dn. (36), 


which is the Weyl correspondence rule. 
The Wigner function is now defined as the Weyl transform of the density 


matrix operator (a factor(1/27f) is added in order to use simple integration over 
the phase space): 


fur(ase) = seq feta — 0/218) + 0/2}e 


or 


fw(a,p) = 5 i e"P/8 6°(q + 0/2) 0(q — 9/2) dn. (37) 


This is equivalent to the definition given in Bq.(8) 
Several properties of fyy echo its origin of classical analogue of the density 
matrix. In particular ; 


/ fw(q,p)dp = |4(4)/?; (38) 


| fw(ap)dg = |e(p)}?, (38) 


where ¢(p) is the wavefunction in p representation. Therefore, if F(q) and G(p) 
ate functions only of q and p, respectively, then 


| [ Feotwtap dade = (F) (40) 


i [ Gtortwr(as2)éadp = (6) (41) 
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More generally, if A is an observable function of q-and p, it can.be shown that 


(A) = Tr(pA) = / / A(q, 2) fw(q,p)dadp, (42) 


where A is the Weyl transform.of A. 

On the other hand, it is not possible to give to fy a simple probabilistic 
interpretation for the reasons indicated at the beginning of this section, and this is 
confirmed by the fact that fw takes in general negative, as well as positive, values 
in phase space. 

A physical insight of the meaning of the Wigner function can be obtained by 
_ observing that it has. large values in the regions of q-space where there is-a strong 
autocorrelation of the wavefunction, and for p values corresponding to the Fourier 
components present in this correlated parts of the wavefunction. Thus, fw has 
large values in regions of phase-space where the presence of the particles can be 
“felt” within the uncertainty principle, and with some more information related 
to the quantum mechanical phases, that can result in negative values for fw. 

Concluding this brief introduction to the Wigner function, we note that if 
the interaction of one particle with the rest of the system can be described by a 
perturbation potential V(q), an equation for the Wigner function can be written 


(Toda et al., 1983), 


Ofwla?, Ofw(apt) _ 1 h 0 Le 
“tela t +2 Tene == (¢- 5235) - V+ ai op twat) (43) 


which reduces to the classical Liouville equation for i + 0. 


A more explicit form of the above equation can also be given, which is more 
suitable for application to real systems: 


Ofw(a,p, t) 4 2 Ofwla pst) : 
at m  @Oq | 


= aay [oP if dye (V(e + by) -V(e—Sullfwlesp +P). (4) 


In order to extend the above theory to the case where scattering is present, an “ad 
hoc” collision term (24¥-)..1: is added to the rhs of Eq.(44) (Iafrate, 1988), but it 
may not necessarily express the same phenomenology as the corresponding term 
in the Boltzmann equation. : 

Some interesting results can be obtained by considering moments of Eq. (43) 
in relaxation-time approximation. Equations are obtained which reduce to the mo- 
ments of the Boltzmann equation in the semiclassical limit, and contain quantum 
corrections which can be included in the classical picture (Iafrate, 1988). 

In what follows we shall report some recent work performed with the WF 


approach on quantum systems of particular interest in modern microelectronic 
research. 
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A. Wigner function applied to the study of resonant-tunneling devices 


Several theoretical works devoted to the development of the Wigner formalism 
for quantum electronic transport have appeared in the literature (Barker, 1980; 
Barker and Lowe,.1981; Barker and Murray, 1983; Lin and Chiu, 1984; Barker, 
1985; Frensley, 1987; renalee. 1988; Kluksdahl et al, 1987, 1989). As a general 
comment, however, it can. be observed that most of the literature on the WF for 
transport problems is primarily concerned with the formulation of the problem, 
rather than its solution. In particular the complicated form of the collision term in 
_ the quantum kinetic equations requires. either severe approximations or formidable 
numerical effords. 

When quantum ballistic systems are considered collisions with phonons/im- 
purities are substantially reduced by a combination of very high mobility materials 
and short channels.(Barker and Murray, 1983). In this case the complexity of the 
collision term can be highly reduced by suitable approximations, and the discussion 
of the quantum features becomes easier. 

The first attempt at calculating the WF for an actual physical system has been 
presented by Barker (1985) and it considers propagation of an incident electron 
gaussian wave packet on a single quantum well, with two bound states EZ, and Eo, 
at the center of a very wide barrier, thus forming two adiacent wide barriers. 

Fig. 7 shows snapshots of the position probability distribution (a), momentum 
distribution (b), and Wigner distribution (c) for a gaussian wavepacket incident at 
resonance when the momentum width of the packet is greater than the resonance 
peak in the transmission coefficient. The complex central structure in the WF 


originates from correlation between the reflected and transmitted wave in phase- 
space. 


The same quantum structure has extensively been studied, using the WF, as 
a resonant-tunneling diode (RTD), (Ravaioli et al., 1985; Kluksdahl et al., 1987; 
Frensley, 1987; Frensley, 1988; Kluksdahl et al., 1989), The quantum well consti- 
tutes a resonant-tunneling system, with a resonant energy marked by preferential 
tunneling. A common step forward of these works is the inclusion of the role of the 
contacts. They are described as “ideal” infinite reservoirs of thermally-distributed 
carriers which act like a source of injected randomly-distributed electrons into the 
device at one terminal, and a perfect sink absorbing all incident carriers random- 
izing their state without reflection upon the other terminal. These contacts serve 
as a boundary for quantum correlations , they remove size dependencies , and 
introduce time-irreversibility in the dynamical evolution of the carrier system. 

Ferry and coworkers (Ravaioli et al., 1985; Kluksdahl et al., 1987; Kluksdahl et 
al., 1989) consider a self-consistent modelling of the RTD based upon the coupled 
solution of the WF equation of motion to the Poisson equation, and include elec- 
tron scattering through a simple relaxation-time approximation. The equations 
are solved with a finite- difference approximation scheme. 

A particular analysis has been devoted to the choice of the initial WF. 

Fig.8 shows the WF for a gaussian wave packet interacting with the resonant 
quantum potential barriers at different times. Quantum interference and tunneling 
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Following this approach the I-V characteristics of the device is obtained by 
increasing the bias potential to its maximum, and then decrementing it ‘towards 
zero, with the current being calculated on the way. The resultant curve shows an 
intrinsic bistability confirmed by the experiments. 

For bias conditions near the peak of the I-V curve an overall depletion of elec- 
trons in the cathode.occurs, as it appears from an analysis of the WF, reported in 
Fig.9-a for this case. As the bias increases further, the resonant. charge becomes 
evident in the WF, illustrated in Fig. 9-b. In Fig.10 the difference between the 
steady-state WF for increasing and decreasing potentials in the bistable region is 
plotted. For decreasing potentials less current flows through the resonant struc- 

- turé, and the injected carriers accumulate on the cathode side of the structure, 
shown as the peak in Fig.10. © 

This analysis of the RTD with the use of the WF has been extended to 
transient regimes (Ravaioli et.al., 1985, Frensley, 1988), the frequency-responce of 
the RTD (Frensley, 1988), and of particular features, like the anomaly in the I-V 
curve for very low applied bias (Kluksdahl et al., 1989). 

The difficulty in extending the formalism to 3-dimensional systems is mainly 
related to computer limitations in memory and CPU time... 


THE PATH-INTEGRAL APPROACH 


There are essentially two different ways of using the path-integral approach 
for the evaluation of the density matrix. One of them is based on the starting 
expression for the evolution operator given in Eq.(14). It is called the “real-time 
path-integral approach”, and we shall discuss it some-more at rengrh later in this 
section, reviewing also one significant application. 

The second way is called “immaginary-time path-integral “approach”, and is 
based on the formal analogy between the analytical form of the equilibrium densi- 
ty-matrix operator, 


-—SH 
? 


po xe (45) 


where 6 = (KT)~?, and the evolution operator for a time-independent hamiltonian 


U(t,0) = e7* Mt ~ (46) 


From the comparison of the last two equations it is easy to understand that there 
can be a formalism (see Feynman and Hibbs, 1965) where the path-integral ap- 
proach to the evolution operator can be extended to the evaluation of the equi- 
librium density matrix with paths along an immaginary time proportional to the 
inverse temperature. Since this method, however, yields only equilibrium proper- 
ties of the system, we shall not discuss it here in details. We simply mention the 
work performed by the North Carolina group (Register et. al., 1988b). They eval- 
uate the equilibrium density matrix for some particular potential profiles related 
to tunnelling problems, adding a stochastic potential to model the effect of the 
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phonons. Their results suggest, as expected, that the presence of phonons may 
reduce the tunneling process by reducing the phase colierence of the-electrons. 

As indicated above, in the real-time path-integral approach the evolution 
operator as given by Eq.(14) is applied to the density matrix. The resulting 
expression 


yt q' st ; ‘ 
Naa a= i da i de J | Dar) i Dal (r eH Sla(-S(@" Ml p(qe, a! 2) 
5 te q 


fits 
(47) 
may be elaborated in a useful way by factorizing the effect of “external agents” 
with respect to the system of interest. In our case we call.z the variables of the 
” “small” system of interest (for example an electron) and X the variables of the 
interacting “large” system (for example the phonon bath). The exponential in 
Bq.(47) can then be factorized as follows: 


ut ek S2)— 52") 0 HIS(X)—S(X)+5(2,X)-S(2',X WDa(r)Dz' (r)DX(r)DX'(r), 


(48) 
where S(z) and S(X) are the actions for the small and the interacting systems, 
respectively, and S(z,X) the action of interaction between the two systems. The 
integrals over phonon paths involve only the second exponential, and an “influence 
functional” can be defined (Feynman and Vernon, 1963), 


Xt pX',t : : Foe 
Farha = fC" fl” Dx (eye (r}\tlsOr-sa4S12-39-Se.491, (49 


such that the evolved density matrix is written as 


p(q,q') = : dq: rt dq; / Da(r) i Da! (r)F(2(7),2'(r))ee@)-S | (50) 


Here for any given path of the system of interest, F carries the information of the 
influence on that path of the integral of all paths of the interacting system. The 
theoretical step performed by introducing the influence functional is not trivial 
since it includes all the effects of the interacting system influencing the behaviour 
of the system of interest. ; 

However, the explicit evaluation of the influence functional is in general pro- 
hibitively difficult, and approximations must be made as in the more standard 
formulations of the problem. For systems where the coupling action is a linear 
function of the coordinates of the interacting part and for systems weakly cou- 


pled, an analytical evaluation of the influence functional is possible (Feynman and | 


Vernon, 1963); for general systems however this is not-true. 

The real-time path-integral technique was applied, without recurring to the in- 
fluence functional, in the pioneer work of Fischetti and DiMaria (1985). The effect 
of the phonons was included by introducing in the action a self energy obtained 
by solving iteratively the Dyson equations for the electron-phonon interaction. 
Electron paths were then sampled with a Monte Carlo algorithm. 
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As amore recent application of the real-time path-integral method to electron 
transport in semiconductors, we will refer to a work (Mason and Hess, 1989) 
that attacks the problem of the transient electron response in a homogeneous 
semiconductor after the application of a constant and uniform electric field. In 
order to be able to handle the influence functional analytically, the authors take 
a linear model for the electron-phonon interaction. In the resulting expression 
another approximation has to be made for the spectral density of the phonons. 
Finally, when the electron paths are sampled with a Monte Carlo technique, the 
space of possible paths had to be reduced to a manageable size, and for this purpose 
the time of integration had to be reduced to the order of one semiclassical collision 


time as for the applications of the QMC technique for the density matrix described 
above. 


The problem seems to be the same in all such quantum-transport numerical 
calculations: multiple integrals with strongly oscillating kernels need to be eval- 
uated. Any naive sampling technique is bound to fail because it samples very 
large values that should cancel with each other leaving a final result many orders 
of magnitude smaller. A smart technique must be found to handle such strongly 
oscillating integrals if numerical results are sought for times at which steady-state 
conditions are reached. Some attention may be deserved, in this connection, by 
the windowing technique proposed by Register et. al (1988). 

An.example of results obtained in (Mason and Hess, 1989) is shown in Fig. 
11. In part (a) of the figure the normalization of the density matrix is shown as 
a function of the collected statistics. It does not reach unity even at the highest 
statistics because of the limitation in the path space sampled with the Monte 
Carlo algorithm. If a larger path space is sampled, the results for average physical 
quantities are in principle more exact but are affected by a bigger statistical error. 
In part (b) of the figure the average space covered by the electrons during the time 
¢ is given again as a function of the amount of collected statistics. From its limiting 
value an average drift velocity (averaged over the past) can be obtained. However, 
owing to the approximations indicated above no comparison can be made with 
equivalent semiclassical results. 


GREEN-FUNCTION APPROACH 


We shall not discuss here in further details the Green-function techniques since 
several speakers will deal with this subject. For numerical applications, in particu- 
lar, we refer to Reggiani’s lectures on “Monte Carlo algorithms for nonequilibrium 
quantum transport”. 

It is worth mentioning, however, the work of the Houston group (Lei and 
Ting, 1985, Xing and Ting, 1987), The old idea of momentum and energy balance 
equations with a Maxwellian distribution is resumed in their work and transferred 
to a Green functions quantum approach. The time evolution of the density matrix 
is written to first order in the electron-impurity and electron-phonon interaction, 
starting from an initial condition that already contains, as parameters, an electron 
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drift and an electron temperature different from the lattice temperature. By appli- 
cation-of Green function tecniques balance equations are obtained whose solution 
yield the electron drift velocity and temperature. 
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Fig.10 : 


Fig.11: 


diagram of a seventh-order. term of-the perturbative expansion for the elec- 
tronic density matrix including electron-phonon and electron-impurity inter- 
actions (see text). 

quantum distribution function as a function of k along the field direction 
containing terms up to the fourth order (continuous line), compared with 
the corresponding classical one (dashed line). The third curve (dotted line) 
represents the ballistic translation of the initial distribution function (Brunetti 
et. al., 1989). 

quantum electronic distribution function as a function of k along the field 
direction containing terms up to the fourth order (continuous line), after 50 
fs from the initial conditions. The dashed curve is the classical distribution 
function at the same perturbative order obtained from the Boltzmann equa- 
tion. The dot-dashed curve is the initial distribution after the ballistic shift 
produced by the electric field (Menziani et. al., 1989). 


: classical electron distribution as a function of energy for a a simplified GaAs 


model at t=100 fs after excitation. The highest peak at 1000 K is the initial 
distribution at t=0 (Brunetti et. al., 1989). 


5: quantum electron distribution as a function of energy for the same simplified 


GaAs model as in Fig.4 at t=100 fs after excitation (Brunetti et. al., 1989). 


: comparison among the CMC drift velocity (continuous curve), its perturbative 


expansion up to the fourth order (dot- dashed line), and the quantum fourth- 
order result (Rossi and Jacoboni, 1989). 


7: momentum (a), and position (b) distributions at a short time after the col- 


lision of a gaussian wave packet with a resonat-tunneling barrier. (c) shows 
the Wigner distribution of the system at the same observation time ( Barker, 
1985). 


: gaussian wave packet interacting with resonant quantum potential barriers. 


The barriers are indicated by the dark band. (a) The incident wave packet, 
moving from left to right, is just beginning to interact with the barriers. (b) 
Gaussian wave packet during reflection. The incident and reflected compo- 
nents are visible, as is the correlation centered around k=0. Part of the packet 
is tunneling through the barrier. (c) Gaussian wave packet after reflection. 
Most of the wave packet has been reflected. The tunneling packet is visible 
to the right of the barriers (Kluksdahl et. al., 1989). 

(a) steady-state Wigner distribution at the peak of the I-V curve: Depletion 
in the cathode region is evident. At the left contact the incoming distribution 
appears as a shifted Fermi-Dirac distribution. (b) steady-state Wigner dis- 
tribution at the valley of the I-V curve. Depletion is strongly evident in the 
cathode region. The distribution in the cathode -barrier-interface region forms 
a quantised state, the ring structure to the left of.the barriers (Kluksdahl et. 
al., 1989). 

difference between the bistable Wigner distributions at bias of 0.36 V. The 
quantised state in the cathode well has more carriers. More current is flowing, 
indicated by the “ridge” in the distribution (Kluksdahl et. al., 1989). 
expectation values as functions of the number of paths sampled. Sampling 
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of each array includes 250.000 paths. (a) shows Tr(p), (b) represents @, af 
is an equivalent scattering rate that gives the high-temperature friction, and 
A is a reference-length: \ = ./h/my. All are fort = e and KpT = 1.2hy7 
(Mason and Hess, 1985). 
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NATO ASI on Granular Nanoelectronics 
Il - Ciocco, Italy 


23. July - O04. August 1990 


RAMAN SPECTROSCOPY FOR THE CHARACTERIZATION 
OF LOW DIMENSIONAL SYSTEMS 


Gerhard Abstreiter 


Walter Schottky Institut 
Techn. Univ. Miinchen 


D-8046 Garching 


—— 


The following aspects are discussed in the lecture 


Growth aspects for the achievements of high quality Si/Ge 


heterostructures and superlattices 


~ critical thickness 

- interface sharpness 

~- phonon Raman spectroscopy as a sensitive local probe of 
various interface and superlattice properties 


~- electronic and optical properties 


Electronic excitations in low dimensional systems with 


emphasis on GaAs based heterostructures 


~- basic concepts of single particle and plasmon excitations in 


3,2,1, and 0 dimensional systems 


q-dependence of electronic excitations in various dimensions 


~ experimental light scattering techniques 


discussion of available experimental results 


Recent advances and new approaches for the achievements of 
lateral microstructures - direct analysis with spectroscopical 
techniques with high spatial resolution (Raman and 


luminescence) 
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Abstract. Experimental results of electronic excitations in microstructured GaAs 
quantum wells are discussed. Micro—Raman spectroscopy and grating coupler effects 
are used to extend the accessible wave vector range of in—plane single particle and 
plasmon excitations in multi-layered two-dimensional electron systems. In lateral 
quantum wire structures a strongly anisotropic dispersion of the plasmon resonances 

is observed, which can be understood in terms of plasmons confined normal to the 

wires. , 


1. Introduction ' 


The first observation of inelastic light scattering by free carriers in semiconductors 
was reported by Mooradian and Wright [1] in doped n-type GaAs more than twenty 
years ago. In this pioneering work plasmons of a degenerate electron gas, coupled to 
LO-phonons were studied. Shortly afterwards Mooradian [2] also reported light 
scattering by single particle excitations in n—GaAs. The ‘early work stimulated 
extensive light scattering research of electron and hole plasmas in semiconductors (for 
reviews see Refs. (3,4,5]). Especially resonant inelastic light scattering experiments [6] 
showed that the method is sensitive enough to observe electronic excitations with 
densities of the order of 10!! cm-2. This was followed by the proposal [7] and the first 
observations of resonant light scattering by quasi-two-dimensional (2D) electron 
systems in selectively doped GaAs—(AlGa)As heterojunctions {8] and multiple 
quantum wells [9] about ten years ago. Numerous publications have appeared since 
then, which demonstrate the versatility of light scattering for the study of various 
properties of 2D carrier systems. These include single-particle and collective 
intersubband transitions in various potential wells which allow a separation of the 
depolarization shift, 2D plasmons in layered electron systems and in-plane single 
particle excitations. The light scattering work of 2D-systems of the past ten years has 
been reviewed for example in [10,11]. So far the dispersion of the different excitations 
could be studied only in a limited wave vector range (q < 1.7 x 105 cm!) in normal 


backscattering geometry from the surfaces. In the present paper we first review 
shortly new developments in electronic light scattering by 2D systems with large 
scattering wave vectors which become accessible by micro—Raman scattering [12] and 
by using grating coupler effects (13,14]. In the last part we concentrate on first light 
scattering results in wire structured electron systems in GaAs quantum wells [15,16]. 


2. Electronic excitations in systems of different dimensionality 


Elementary excitations of free carrier systems fall into two main categories: collective 
excitations and single particle excitations. Resonant electronic light scattering allows 
a separation of the two kinds of excitations by analyzing the polarization of incident 
and scattered light. Polarized spectra (incident and scattered light polarized parallel 
to each other) usually display collective excitations, i.e. plasmons which are coupled 
to LO—phonons in polar semiconductors. Single—-particle excitations are observed in 
crossed polarizations. They are characterized by uncorrelated excitations of electrons 
below the Fermi energy to empty states above the Fermi energy. In backscattering 
geometry the scattering wave vector is given by q~2-2an/A, which is about 
7-105cm7! for GaAs in the red spectral range (7 is the refractive index and A the 
laser wavelength). It leads to the kind of excitations shown in Fig. la for an isotropic 
3D-electron system. The cutoff frequency of the single particle excitations is at about 
G°Vq where Ve is the Fermi velocity. In 3D they have approximately triangular 
lineshape. The plasmon frequency is proportional to Vv O35 , with n,» being the 
three—dimensional carrier concentration. It exhibits a weak dispersion which in 
random phase approximation (RPA) is given by w2(q) = w2(OS + 3q2 ve/5. 

In 2D-systems, which are realized for example in narrow quantum wells or 
heterostructures, the motion of the carriers is quantized normal to the layers and is 
free in parallel direction (q,). Typical excitations of such systems are shown in 
Fig. 1b. In usual backscattering geometry normal to the layers (q,= 0) collective and 
single particle intersubband excitations are observed, which reveal the subband 
splittings wo, and the depolarization shift. Recent experiments on extremely high 
mobility samples seem to indicate that depolarized spectra (scattering via spin 
density fluctuation) do not give exactly the one particle energies wo1, but are subject 
to measureable exciton—like shifts to lower energies [17], which were neglected in most 
of the earlier experiments in GaAs. In—plane single~particle and collective excitations 
are only observed for finite q,. Due to the high refractive index of GaAs the 
maximum achievable in-plane wave vector in backscattering from the surface is given 
by q,, = (47/A) sina < 1.7 x 105 cm‘ for laser energies around the Eo and Ep + Ao 
energy gap. The angle a is defined between the incident (scattered) light and the 
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Fig. 1 Single particle and collective excitations in 3, 2, 1 and 0 dimensional 
electron systems (schematically) 


normal to the sample surface. The dispersion of both plasmons and single particle 
excitations of layered electron systems have been studied in this limited wave vector 
range [10,11,18—21]. An extension of the accessible in-plane wave vector is discussed 
below (see also [12,13,14}). 

In Fig. ic and 1d the possible excitations in 1D and OD systems are shown 
schematically. In quantum dots, confinement exists in all three space directions which 
finally results in an atomic—like situation with no dispersion in any direction. In 1D 
systems (quantum wires) free motion of carriers is still possible along one direction 
(qy). For the other two directions different confining potentials are assumed in 


Fig. lc. This leads to different single particle and collective excitations in all three 
directions. The restriction of free motion in one dimension also changes the single 
particle excitation along the wire drastically. No low frequency excitations are 
possible for a finite wave vector qy along the wire. First experimental results of 
electronic excitations in wire structured quantum wells or heterostructures have been 
published recently [15,16]. Our results, especially on confined plasmons, will be 
discussed in this paper. To the best of our knowledge, no results have been published 
so far for quantum dots. 


3. In—plane excitations in layered electron systems 


Until recently in—plane excitations in 2D systems were observed only in a relatively 
small wave vector range. In backscattering from the surface of a layered electron gas 
the total wave vector is given by q = 2-2an/A, which divides into q, = (47/A).sina 
and q, = (4xn/A) J 1-sin2a/72, depending on the angle a. A maximum in-plane 
wave vector is achieved in back—scattering geometry from the edge of the layered 
structure (cleavage plane). Such geometries have been realized recently by using low 
ter yerature micro—Raman scattering [12]. Results of single-particle excitations and 
plasmon excitations in cleaved and wedged GaAs multilayer structures are shown in 
Fig. 2 and 3. The maximal achievable in-plane wave vector is 
q, = 4an/A = 7-105 cm! using a red laser line for excitation. The sample consists of 
25 GaAs quantum wells with a thickness of 60 A and a carrier concentration of about 
9 x 10!! cm*2 per well. The in—plane single particle excitations shown in Fig. 2 exhibit 
the characteristic lineshapes of 2D systems. The sala, ines are calculated lineshapes 
using the Lindhard—Mermin dielectric function of the 2D electror. systems. All spectra 
can be fitted by assuming a temperature of the electron gas of T po lsk. Scattering 
from the as grown surface (q, = 1.3-105 cm" for a = 45°) and from the cleavage 
plane reveal a carrier concentration of about 9 x 10!! cm-? and an electron damping of 
rs ~1meV. The polished surfaces show somewhat lower quality which results in a 
slightly reduced carrier concentration and larger damping. 

Plasmons in layered electron systems are characterized by the Coulomb 
interaction between the layers. The dispersion of a purely 2D system is given by 
w3(a,) = 2m, ¢%q,/m*e, where Nyy is the 2D carrier concentration and ¢ the 
dielectric constant. The coupling between the layers leads to a plasmon band which is 
determined for an infinite system by u(q,,49) = Wor (4,) 5 (Gyrdg)s Where S (q,,,49) = 
vsinhq, d/(coshq, d—cosqgd) is a structure factor. For a finite number N of electron 
gas layers N plasmon branches are obtained [22] which can be approximated by the 
above-mentioned formula with qg = (21/Nd)8 with 8 = 0,1,...N-1 for a large number 
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Fig. 2 | Depolarized Raman spectra of layered electron systems. The dotted lines are 
theoretical fits to the observed single-particle excitations using the 
parameter given in the figure. The in—plane wave vector is determined by 
the scattering geometries which can be realized with micro—-Raman 
spectroscopy as shown schematically. 

Fig.3 Plasmon excitations (polarized spectra) of layered electron systems 
measured with micro—Raman spectroscopy 

of layers. The branches of the plasmon band which are probed by inelastic light 

scattering depend on the ratio q, /q,- In backscattering from the as—grown surface q, 

is large so the lower (acoustic) plasmon modes are observed. From the cleavage plane 

q, = 9, therefore the more 3D upper branch of the plasmon (B = 0) is probed. We see 

again from the spectra of Fig. 3 that the edges polished under an angle have less 

quality compared to the grown or cleaved surfaces. The plasmon dispersion, however, 
can be measured up toq, = 7 * 105 cm"l. 

Another way to extend the accessible wave vector range is the so—called grating 
coupler effect (13,14). A periodic structure on top of the layered electron system can 
couple larger in-plane wave vector into the underlying semiconductors due to 
diffraction. The total wavevector is then given by q®=q,+m-g_ where 
m=m, + Mg is the diffraction order of incident and scattered light and g = 27/a is 
the reciprocal lattice constant of the periodic grating. Examples of collective in-plane 
excitations in GaAs quantum wells are shown in Fig. 4 [14]. A periodic etch mask of 
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Fig. 4: Plasmon excitations 
measured with resonance Raman 
spectroscopy in layered quasi 
2D electron systems. The different 
plasmon peaks correspond to 
different in—plane wave vectors qy 
which are accessible via diffraction 
of light due to the grating coupler 
which is etched into the sample. 
The numbers label the grating 
order. The 0 order peaks are 
similar to the plasmons observed in 
the upper three spectra of Fig. 3. 
The seventh order corresponds to 
an in-plane wave vector of 
q, = 5.5 x 105 cm"! (from (14). 
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photo resist stripes with a period of 8100 A was produced by holographic lithography. 
The grating was etched into the sample in an optimized reactive ion etching process 
using SiCl, plasma (23). Eighteen quantum wells remained underneath of this grating 
with a carrier concentration of about 5 - 10!! cm-2. Plasmons can be excited up to the 
seventh order now, even in backscattering with normal incidence on the surface. The 
peaks in Fig. 4 are labeled by the suitable diffraction order index m. For a= 0° the 
scattering wavevector is qa= mg. With decreasing a we have to add q, = (41/A)sine 
to obtain the total in-plane wavevector. With this grating coupler technique, the 
dispersion of layered plasmons could be measured up to about q,=5.5 = 105 cm71(14]. 
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Typical geometry of the wire structured samples achieved by holographic 
lithography and reactive ion etching. Each wire contains five modulation 
doped GaAs quantum wells. Details of a wire are shown schematically on 
the right—hand side of the figure. 


4. Electronic excitations in wire structured quantum wells 


In the last section we discuss very recent results on inelastic light scattering by free 
carriers in lateral quantum wire superlattices. The samples are prepared in a similar 
way as the etched grating coupler. The resulting periodic wires are shown 
schematically in Fig.5. The geometrical width of the individual wires is about 
4000 A. Each wire contains five sheets of electrons in GaAs quantum wells with a 
thickness of 60 A and a 2D carrier concentration of about 5 10! cm-2. The actual 
lateral width of the electron systems is expected to be considerably smaller than the 
geometrical width of the wires due to surface depletion. This is shown schematically 
on the right-hand side of Fig. 5. The Fermi level pinning at the surface due to surface 
states and damage causes a depletion layer on each side, which leads to an 
approximately parabolic potential well toward the surface with a flat part in the 
middle [24]. A subband splitting due to this confinement in x-direction of about 1 to 
2 meV is expected. 

We have investigated both collective and single-particle excitations in such wire 
structured system by resonant inelastic light scattering. The dispersion of these 
excitations were studied in backscattering by tilting the laser parallel and 
perpendicular to the wires. In Fig. 6 polarized spectra are shown for three tilt angles 
in each direction. Several plasmon excitations are observed in each spectrum. No shift 
of the peak position is observed with tilt angle a normal to the wires (7 = 0°). On the 
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Fig. 6: Typical plasmon excitations 
in the wire structured sample for 
various tilt angles a normal to the 
wires (y = 0°) and along the wires 
(y= 90°). The peak positions 
remain unchanged for y= 0° (no 
dispersion). With +y=90° a 
considerable dispersion of the 
plasmon modes is observed, which 
ig strongest for the lower order 
modes. 


RAMAN INTENSITY 


ENERGY (meV) 


other hand the peaks shift to higher energies when a is increased in the direction 
along the wires. The dispersion is stronger for the low energy excitations. The 
vanishing dispersion perpendicular to the wires shows that Coulomb coupling between 
neighbouring wires is very weak and we observe resonances determined by the isolated 
wires, The spectrum of plasmon modes obtained parallel to the wires reflects 
qualitatively all the predictions made by Eliasson et al. (25] for a laterally bounded 
single layer 2D electron gas. To include the multilayer structure, we use the simplified 
model of confined plasmons. The interpretation of the plasmon resonances is based on 
the dispersion shown in Fig 7. The dotted lines are the five plasmon branches 
calculated for Don = 5.1 - 10% cm-2, We keep this plasmon dispersion of the layered 
2D electron gas and impose standing wave conditions for the plasmons in direction 
perpendicular to the wires. Then the in-plane component of the plasmon 
perpendicular to the wires is restricted to values qm = m-2/We, where We is the 
effective wire width which takes the partial lateral depletion into account. The 
component parallel to the wires remains unrestricted: q, = ¥ q2 + ay(a)®. These 
values are represented as vertical lines in Fig. 7. The peaks are assigned to the 
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Fig. 7: Dipersion of plasmon 
excitations in the wire structured 
sample for y=0° and y= 90°. 
The allowed values for the 
in—plane momentum for the 
different angles a are represented 
by the vertical lines. The effective 
wire width determined from this 
analysis is We = 2400 A. 
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quantized momenta, for which their energy coincides with one of the plasmon 
dispersion branches. The peaks in Fig. 6 are labeled according to the quantization 
index m. Different branches of the plasmon dispersion are distinguished by primes. 
From the lowest confined plasmon, labeled by 0, we~ obtain a perpendicular 
momentum of 1.3 x 105 cm"! which results in an effective wire width of we x 2400 A. 
These experiments represent the first experimental study of plasmon dispersion and 
plasmon confinement in lateral quantum wire superlattices by resonant electronic 
light scattering. A detailed analysis of the results is presented in [16]. 

In conclusion, we have shown that resonant electronic light scattering is a 
versatile tool to study the dispersion of single particle and plasmon excitations in 
electron gases of different dimensionality in a wide wavevector range. Especially the 
possibilities opened by micro~Raman spectroscopy and grating couplers, and the new 
properties of quantum wires and quantum dots, which become now accessible, will 
lead to stimulating and exciting developments in the near future. 
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Nanofabrication of Quantum Wires and Rings 
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Outline: 
Dimensional requirements 
Lithography for nanostructure definition 
Pattern transfer techniques 


(Direct fabrication - wait for next lecture) 


Issues in Fabrication of Quantum Confinement Devices 


¢ Dimensional! requirements 


+ Preservation or improvement of structural resolution during pattern definition 
and transfer 


¢'Cleanliness' 


- start with structure with high degree of perfection, usually grown by 
MBE. Critical, active regions are buried below surfaces, interfaces are 
of high (~monolayer smooth) quality 


- we do something nasty to this beautifully perfect structure using 
electrons, ions, dopants, etchants, contacts, chemicals....... 


- and then hope that we have a perfect system in which to demonstrate 
some new physics or make smart devices 


Fabrication of Structures for Quantum Confinement 


Starting materials: 
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ONG 


Quantum well 


1. Quantum well layer 


+V (2-3V) 
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OPDGEPATEPTOOOOTEEOPEERPEOLELEDEOBEDET. 


LECLELPIOTIILEELEITEDEDETIITALEELOLELE 

LEILOLPADLODIOTLELIGPIPLLILIGLEOEEOES. 

CULECAILTOTTIOLLESISILPPLELEGLLOLOEEL, 
GELOCLELEETOBIELEBLOTLEPOTEDEEEEEEIOEE, 


c. Deep etching 


a. Gating b. Shallow etching 


The Device Process Cycle 


+ Structure fabrication usually requires a multi-step approach 


Growth 


Lithography 


Direct 


route Deposit 


Passivate 
| 


Package 


* All steps in the process require to preserve or enhance the resolution of their predecessors 
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Lithography for <100nm dimensions 


Odie 4 2DEG quantum 
wires 


Contacts, 
interconnects 
vertical devices 
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devices 


Suitable techniques are 
* x-ray lithography: but this is a contact printing process relying on other 
methods to make masks 


* Ton beam lithography 
* Electron bea lithography 


Photolithography 


Photolithography 
¢ Limited to ~ 0.4 - 0.5um (arguably) 
* But specialist techniques or tricks can achieve 


smaller dimensions e.g. 


- Shadowing 
Angled evaporation 
Resist 
Substrate 
- Holography 


Overlapping laser beams generates 


interference fringes with period 


ho 


Ne 2 sin 8 


| Sal 


aie a 


X-Ray Lithography I 


If diffraction effects are a limitation, move to shorter wavelengths e.g. 
X-rays : 


Problem: no lenses. Must use contact or proximity printing process 


* need some other non-optical method of mask making e.g. E-Beam 


Choice of wavelength: Depends on degree of contact and mask 


transparency 
0.1nm 10nm 
—__———_—_——__ ——— > _ increasing 
Mask (& dust) <————__» _ Mask more opaque 
transparent 
Less diffraction <——_—_—————__» __ Diffraction worse 


Increasing photoelectron 9 <———_-___ > Decreasing photoelectron 
range range 


X-ray Lithography II 


Mask Technology 


Si + Polyimide + Gold Absorber 


Au absorber pattern 
Polyimide 


B-doped Si 2um thick 


Quartz support ring 
Polyimide + gold (conformable mask for contact) 


- Transparent to | Au absorber patiern 
long 


wavelength 
radiation 

¢ Flexible, 
electrostatic 
chucking 


Polyimide 


Quartz support ring 
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Inorganic Resist Techniques for <10nm Lithography 


Attempts at higher resolution resist processes 
Contamination 
¢ Decomposition of vacuum oil or other carbonaceous material 
¢ Carbon deposit can be used as an etch mask 
* 8nm resolution achieved 
Drawbacks 
¢ Surface migration of undecomposed contaminant causes depletion effect limiting 
packing density 
* Known that activated species migrate after interaction with beam: linewidth 
resolution limitation 


Primary beam . 
Contaminant source | Primary beam 


i J fa \iciivated species 


Henlenonst é Migration of activated species 
epletion of contaminant source before ‘fixing’ 


NB Watch out for these effects in proposals to carry out beam-controlled growth or 
deposition. Presumably can control with substrate temperature 
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Inorganic Resist Techniaues for <10nm Lithography II 


Another promising line of attach involves “hole-burning" in inorganic resists. 
Thin films are prepared by evaporation and exposed to very high doses of electrons 
Holes are “burned” in the films 


¢ Examples of materials include magnesium oxide, 
aluminium/strontium/cadmium fluoride (last two are epitaxial on GaAs) 


¢ Resolution <~2nm has been achieved 


* Require very heavy dose (10's of C/em42) for exposure (cf PMMA ~104-4 - 
104-3 C/em‘2) 


* DOSE RATE DEPENDENT (unlike ordinary resists): will not expose at low 
current 


Problem: material does not clear completely in smallest structures 


Practical resolution (after transfer) no better than best polymers 


Interesting nonetheless: basis for all-vacuum process 


Electron Beam Nanolithography 


Electron gun Resolution determined by: 
ewe: (25-100K ¥) + Probe diameter 
(can be 0.5nm) 
” Ys - Electron scattering 
C1 (magnetic lens) (minimise with high energy beam) 


GQ BR » Secondary electron generation 
(38nm?) 
+ Long-range chemistry, 


molecule size effects 


Blanking 


Advantages: 

- Direct computer 
control of pattern 
* High resolution 


Objective 


Scan coils 


Control 
computer 
Digital pattern 
generator 
Drawbacks: 
+ Low throughput 


Substrate on - Limited field size 
laser controlled Let 
stage Applications: 

* Mask making 

+ Direct write <0.5um 


Focussed Ion Beams 


Why Focussed Ion Beams? 


¢ Much less primary elastic scattering of ions than electrons 
especially in resist (less deviation/lateral scattering) 


¢ Can be arrested by resist - no penetration or damage to the 
substrate 


¢ A focussed ion beam offers other fabrication opportunities: 
- Direct Implantation 
- Direct Etching/Sputtering 
- Localised damage for quantum wire isolation/definition 
But: 


* The technology is immature: viable. bright sources have only 
recently emerged 


¢ The resolution is limited by chromatic aberrations (energy spread 
on the beam) to 10-30nm 
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* Sample bombarded by ions from rf excited plasma accelerated across dark 
space by dc bias (20-1000V) 


+ Gas may be inert (sputtering) or chemically active 
+ Chemical reactions enhanced by ion bombardment 
* Products of chemical reactions should be volatile at substrate temperature 


* Resolution = resolution of mask: but damage to material results 


SiCls 


RU 


/ pas 


CisF CCloF2 

He 
$8 only as mixture 
Gas alone 


Sources of Dry Etch Damage 


© Mask penetration 


Desorption 


Bombardment by gas molecules Gidswae daniaue tay 


accumulate 
Effects of direct 
bombardment: 
Disruption Surface compound 
Desorption formation 
Ad/Absorption 
Implantation 


Damage region continuously removed 
on directly etched surface 


Electrical Study of Sidewall Damage 


Damaged region 


10on 


Surface deplet 


Measure conductance versus width 


¢ Perfect wire cuts off at 2 x surface depletion depth 


« Damaged wire has larger cutoff 


ion of the damage depth allows better 
modeling/design of quantum wire structures 
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¢ Character 


impurity Induced Disorder 


impurity into 
quantum well structure 
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2, Implant ions (Ga, Si, B, F) and anneal 


- laser mirrors 


- lateral confinement 
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- quantum wires and dots 


d into GaAs well 
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- Avoid free carrier absorption by use 
disorder 


Damage Confinement Technique 
Method due to Scherer and Roukes, developed by Thornton 
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¢ 2DEG starting material 


¢ Deposit mask (resist, lifted off 
inorganic, e.g. Strontium 
Fluoride 


¢ Bombard with low energy ions 
of noble gas (Ar, Ne) 
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SERSSSSOSSSOSOO SOCIO * Confinement by mobility 
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¢ If carefully controlled, 
lithographic width = electrical 
width 


Fabrication using the Scanning Tunneling Micrascope 


Tunneling current ees 
between tip & 8 
substrate oe 


Tunneling current can be used to: 
« Dissociate (expose) resist 
* Deposit contamination 
- Dissociate metallorganics 
* Directly modify surface 
Possibility of atomic scale fabrication?? 


Conclusions 


¢ Quantum confinement structures: straddle a size regime where 
fabrication is difficult 


¢ Cleanliness, or structural perfection, is a major issue 


- In doped structures, depletion from surfaces and interfaces 
may help to confine electrons from imperfect regions. 
Transport devices can be very ‘clean’ (eg. gated point 
contacts) 


- Depletion can also be utilised to obtain an extra degree of 
confinement: the drawback being the difficulty of exercising 
fine control of the shape and size of the active region due to 
fringing effects 


- In optical quantum structures in which carriers are 
photoexcited, there is greater sensitivity to processing effects 


¢ ‘Traditional’ processing still works quite well for transport 
experiments: more ingenious techniques may be needed for optical 
devices 


Optical Properties of Arrays of MicroDevices 
(Quantum Wires and Dots) 
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Outline: 
Dimensional requirements 


Studies of etched and overgrown quantum wires & 
dots: the role of surfaces 


Evidence for quantum confinement 
Clean fabrication techniques 
Ingenious experiments 
Conclusions 


Introductory 


Increasing confinement from 2- to 1- to 0-D 
Typically by starting from 2-D structure such as QW 
Pattern laterally to create confining potentials: 

- surfaces (etch) 

- modify alloy composition (ID) 

- apply external potentials 


S(E)} S(E) 
Optical transitions restricted 
to few discrete energies 
- enhanced oscillator 


strength 
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Modifications to DOS 


Size constraints 


(After Vahala, J. Quant. Electr. 24,523,1988) 
Confinement 


Lower limit : Energy 


If D too small, no confined sea in dot 
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Upper limit 
If D too large, transitions smeared by thermal 


excitation 
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Also detectability and fabrication tolerance limitations 
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( Conclusiong- Free Standing Kicked Dots 


+ Majority of data points scattered about relalive iniensky of 1 
+ Luminescence scales approximately with volume of quantum 
well material excited 
- No blue shift - dots not small enough 
e INDEPEMVENT OF PROESSING GNRTUNS (MASK ETCH Rivre.) 
issues 
« What is role of surface recombination? 
¢ Unimportant here? 
¢ Dots etched with Ar show low luminescence efficiency 
- Chemical passivation effect? 


- Linear incident power dependence supports low surface 
recombination rate CZ Dimaisves LE FoR subspace LécomBatTiow 


—» + Can carriers diffuse to surface? ~<— 
- Variability 


« Mask variations - cathodoluminescence 
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FIG. |. PLE spectra at zero magnetic field of a QWW with 
t= 70 nm and of the corresponding reference QW. The spec- 
(ra are shifted vertically with respect to cach other for clarity 
and show the dependence on the polarization of the exciting 
laser light. For s-polarized (p-polarized) light E is parallel 
(perpendicular) to the wires. The excitation intensity was 
about 100 mW/cm?. 
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FIG. 2. Magnetic-field-dependent PLE spectra of a ¢ = 7( 
nm QWW. (b), (c), and (e) denote the hh, hhiz, and ty 
transitions, respectively; the index / labels the order of the 
inter-Landav-level transitions. The spectra are shifted verti- 
cally with respect to each other for clarity. The excitation im- 
tensity was about 100 mW/cm?. 


FIG 1 Schematic cross-sectional drawing of the quantum wire structure. 
Parameters used in the calculation of confined energy levels and wave func- 
tions are indicated. 
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FIG. 2. Low-temperature PL (dashed hnes) and PLE (sold hnes) spectra 
of ~ 350-A-wide quantum wires. 
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3. Emission spectra of quantum wire diode lasers with different Al 
fraction yin the cladding layers and different guiding layer thickness ¢ 
Fig. 1). (a) y = 0.5, ¢ = 0.2 um; lower: 40 mA, upper: 50 mA (cavity 
hL=2.07mm,/, = 41 mA). Lasers with longer cavities (as long as 
im) lased at same wavelength. (b) y = 0.7, ¢ = 0.15 um; lower: 22 mA, 
1125 mA (L = 3.48 mm, /,, = 23 mA). 
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FIG 4 Evolution of emission spectra of a quantum wire semiconductor 
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a) Schematic sketch of the quantum well sample and wire mask. The typical wire widths 
the order of 100 nm. b) Lateral distribution of defects due to ion straggling under the 
1 the quantum well plane. 
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‘FIG. 1. Schematic not to scale of the structure after implantation and an- 


nealing. The mask position during implantation is also indicated. The wires i 
are actually ? jem apart. 
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fog FIG 2 Cathodoluminescence spectra (7 < 15K) with the electron probe 
oT? fer? Gee LOPANTATION ( ~ 200 A in diameter) on: a large masked area [(a) left], a large interdif- 


fused area [(2) right}, quantum wel wires with mask sizes of 4500 A (b). 
¢ Gwe Smer 1700 A (c), and 1400 A (d). (€) is a CL micrograph obtained by selecting 
the luminescence energy (3.647 eV) of the QW wires which appear as 


« NEU FRATORES bright areas. The mask size for these QW wires is 600 A. 
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FIG, |. Schematic diagram of the stressed wire and substrate, conduction- 
and valence-band modulation, and exaggeration of the red shift in the band 
gap under the wire. 
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FIG, 2. Phototuminescence spectra of the InGaAs quantum well underly- 
ing (a) an wnetched portion of the sample, (b) a uniformly etched region. 
and arrays of (c) 400-nm and (d) 200-nm-wide wires, and (¢) 400-nm- 
wide (square) dots. All regions of the sample shown here were located with- 
in a few hundred microns of each other. Relative intensity scales are shown 
at the bottom left of each spectrum. Note the ~9 meV red shift of the 400 
nm wires, and slightly smaller and larger shifts for the 200 nm wires and 400 
nm dots, respectively. Sharp spikes on the high-sensitivity spectra are from 
cosmic rays, to which the Ge detector is sensitive. 
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Biosynthesis of cadmium sulphide Zinc-blende ee 
quantum semeconductor 
crystaBtes 


C. T. Dameron*, 8. N. Reese", R. K. Mehra*, 
A. R. Kortant, P. J. Carrot, M. L. Steigerwaldt, 
L. £. Brust & D. R. Winge* 


Counts s** 


* University of Utah Medical Center, Salt Lake City, Utah 84132, USA 
t AT&T Bell Laboratories, Murray Hill, New Jersey, 07974, USA 


NANOMETRE-SCALE semiconductor quantum crystallites exhibit 
size-dependent and discrete excited electronic states which occur 
at energies higher than the band gap of the corresponding bulk 
solid’, These crystallites are too small to have continuous energy 
bands, even though a bulk crystal structure Is present. The onset 
of swch quantum properties seis a fundamental limit to device 
miniaturization in microelectronics®, Structures with either one, val 
two or ail three dimensions ca the nanometer scale are of particular 0 20 a0 60 Bo 
interest in solid state physics’. We report here our discovery of a 
the biosynthesis of quantum crystallites in yeasts Candida glabreta = FG.1 1.54 A powder X-ray patterns of desalted lyophilized Cd- y-Glu pept. 
and Schizosaccharomyces pomode, cultured in the preseace of cad- complexes from C. glabrata (top) and S. pombe (bottom). The samples we 
miu sels. Short chelating peptides of general structure (y-Glu- Purified bya combinatlon of lon-exchange and gel-filtration chromatograr 
Cys),-Ghy control the mactention und growth of CaS crystallites as described previously : . The samples were desalted by chromatograp 
to peptihe-copped ietraceddalar particies of diameter 20 A. Treese on Sephadex G-25 equilibrated and eluted with delonized water. Desait 
quantum CS crystallites are mere ‘i than CdS par- samples were dried in X-ray tubes under vacuum. Solid and dashed tint 
. i 7 expected pattems for zinc-blende and rock-salt lattices, respectively, wi 
ticles synthesized chensicaily. X-ray data imdicate that, at this ong A coherence length. The sha‘p peak at low angle (4°) represents partic 
small size, the CdS structure differs frem that of bulk CdS and 


spactng in the dense Isolate, The nearest-neighbour distances are 22 ha 
tends towards a six-coordinate rock-sek sxrectnre. 20 A for isolates from C. glabrata and S, pombe respectively. 
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ASI - Granular Nanoelectronics 
MOLECULAR ELECTRONICS 


J. P. Launay, Molecular Electronics Group, 
CEMES / LOE, 29 rue Jeanne Marvig 
31055 TOULOUSE CEDEX France 


Abstract. A review is presented of "Molecular Devices" which have been already 
realized. This includes molecular wires, rectifiers, photodiodes, switches and 
bistable elements. Then some basic problems encountered in Molecular 
Electronics are evoked, in particular the need to associate a large number of 
components, the energy problem, and the place of Molecular Electronics in 
information processing systems. Several important research themes are listed, 
including nonlinear behaviours, logical functions at the molecular scale, the need 
to cascade components and to investigate the behaviour of a Quantum Computer. 
Finally some possible short term realizations are indicated. 


INTRODUCTION 


Molecular Electronics can be broadly defined as the processing of electric, magnetic, or optic 
signals with molecule based devices. After the discovery of conducting polyniers in the 70's, 
the idea that molecules could be used to process signals was introduced by Aviram and Ratner 
(Aviram, 1974), who introduced the molecular rectifier concept, and later by the (late) F.L. 
Carter. From the beginming of the 80's, in a series of fascinating multidisciplinar workshops 
(Carter, 1982, 1987, 1988), he launched definitively the concept of Molecular Electronic 
Devices (MED). Of course a number of ideas evoked there were (and still are) highly 
speculative. But this had the immense advantage to wigger an increasing activity in chemical 
synthesis. There are now several example of "Molecular Devices” emerging from the 
laboratories. Most of them are based on simple and often naive analogies with actual Electronic 
Devices. But their fundamental study will certainly be highly instructive and full of surprises. 
The question now is to try to foresee the behaviour of large assemblies of such devices. This 
problem could have strong relations with the concept of Granular Nanoelectronics and both 
approaches (one from the above, the other from the bottom) are destinated to meet somewhere. 


The present paper is organized as follows. In (I), we present the state of the art about 
"Molecular Devices” which are classified according to simple concepts coming from electrical 
analogy. They constitute thus the “spare part box" for future engineers. In part (II) we address 
the basic problems encountered in Molecular Electronics, which all come frome the extreme size 
reduction. Part (Ill) will suggest several research tracks which (we believe) could be crucial to 
transform Molecular Devices from simple Laboratory curiosities to actual useful systems. Most 
of these tracks involve long term research ; however in part (IV) it will be seen that short term 


ee eee te, 


realizations and applications are nevertheless possible and will be extremely useful to reinforce 
credibility for this topic 


1.A REVIEW OF EXISTING "MOLECULAR DEVICES" 


A large number of molecular electronic devices have been suggested in the literature. Here we 
shall restrict to systems which have been actually synthesized, or which appear around the 
corner from a chemical point of view. We furthermore restrict on molecules based on electron 
transfer. Composents based on the propagation of other particles are conceivable but either the 
chemistry is extremely difficult (case of ions for instance) or the transfer process is not 
conservative (case of energy transfer for instance). In addition electron transfer is a universal 
process, already extensively used in conventional electronic devices. It can in principle survive 
to extreme size reduction and there is no theoretical impossibility to perform single molecule 
addressing by electrical means. 


LL_ Wire 


This is the most simple element the fuction of which is to transfer a bit of information 
(represented by an electron). However at the molecular scale it is not possible to study 
independently the behaviour of the different parts of a quantum system. Thus the "molecular 
wire” must encompass not only the "conducting" part, but also the terminal sites the function of 
which is interfacing with the outside world or with another component. 


A first realization makes use of a binuclear coordination complex M-L-M (M = metal site, L = 
bridging ligand) in which the two M sites are able to exist in two oxidation states (e. g. Ru 


(NH3)5 2+/3+) and L is a x conjugated system (Woitellier et al 1989). When the complex is in 
the "mixed valence” state, i.e. one of the metal site is in oxidation state 2+, and the other in 
oxidation state 3+, there is a possibility of electron transfer between the two remote sites. This 
gives rise to a special electronic transition, the intervalence band, which usually occurs in the 
near IR. This band gives some information about the degree of through bond coupling between 
the two sites. Thus one has indirectly some information about the rate of the thermal process : 


Ru2t- L- Rust - Ru3t-L- Ru2t 


Figure 1. A model of rnolecular wire 


It should be noted that this process occurs more or less rapidly through any type of organic 
ligand. At this scale, there is no perfect insulator and the electron transfer occurs through 
tunneling without an actual oxidation or reduction of the ligand. To use semi-conductor 
language, Ru 2+/3+ sites can be considered as impurities for which the energy levels fall in the 
gap between the valence band and the conduction band. The electron transfer is due to an 
indirect coupling between these Icvels, without creation of charge carriers in the conduction or 


valence bands. Recent results obtained in our group show that for the family of complexes 


bridged by a&,@ bipyridylpolyenes, the electronic coupling decreases very slowly with distance 
(Joachim et al, 1990). The longest systems of this type have been described by J.M. Lehn and 
coworkers, for instance a “caroviologen" molecule in which 2 terminal pyridinium sites are 
linked by a chain of 11 conjugated double bonds (Arrhenius et al, 1986). The total length is ca 
34 A and this molecule can span a vesicle membrane and perform electron transfer between the 
outside and the inside (Blanchard-Desce, 1989). Recently, Effenberger has described an 
analogous system but dedicated to energy transfer (Effenberger et al, 1988). 


L2. Rectifi 


In the previous case, the two terminal sites were identical. If now one introduces some 
asymmetry by associating a donor group (D) on one end and an acceptor group (A) on the other 
end, the electron transfer will be favored in one way and rectification properties will appear. It 
should be noted however that the term "molecular rectifier" or "molecular diode" is frequently 
used to designate two different types of molecules. 


The first type corresponds to the original concept by Aviram and Ratner of a molecule with 


topology D-o-A, where o is a bridge containing only o bonds, in order to limit the mixing of 
electronic levels of A and D. This molecule is destinated to be sandwiched between two metallic 
conductors so that it could be acdressed by a permanent electron flux. Several molecules of this 
type have been prepared (Metzger et al, 1986) They were designed for incorporation in 
Langmuir Blodgett (LB) films with the objective of measuring the electrical resistance through a 
single monolayer placed between a Sn O2 and a mercury electrode. This experimental set up 
suffers from two drawbacks : the frequent non perpendicular disposition of molecules in the 
film and the existence of microscopic holes in the films. Some attempts have been made to 
detect rectification with the tip of a Scanning Tunneling Microscope (STM) but they are 
unconclusive at the present time (Aviram et al 1988, 1989). Finally very recently, Sambles and 
coworkers have claimed that a LB monolayer of a potential molecular rectifier sandwiched 
between Pt and Mg electrodes exhibited rectification (Geddes et al, 1990). This seems to be the 
first successful report. 


Figure 2. Examples of molecular rectifiers 


The second type of molecular rectifier uses a m link between D and A sites. This provides a 


strong hyperpolarizability (B factor) and such molecules are currently used to devise organic 
materials for non linear optics (second harmonic generation). In these molecules, contrarly to 
the previous case, the electron transfer is virtual because the D*t-A~ configuration is never 
reached. The effect of the electric field is to alter the degree of mixing between the preponderant 
D-A configuration and the Dt-A~ one (Chemla and Zyss, 1987). 


1.3, Photodiode, 

Molecular photodiodes have been extensively studied in the last year with the prospect of 
achieving artificial photosynthesis. The general principe is to start from a photosensitizer (PS) 
linked to donor and acceptor groups. After excitation, the species PS * has an electron in a high 
energy orbital and a hole in a low lying orbitai (cf electron hole creation in conduction and 
valence bands of a semi conductor). By a proper design of the molecule one can expect an 


electron transfer towards A and a hole transfer towards D, sc that a vectorial charge separation 
is achieved : 


hv 
Depseks. tetas: > D-PS*-A 
Dt- PS - A’ 


Thus with PS = porphyrin, A = quinone and D = carotene (see Figure 3), the exited state has a 
lifetime of 3 psec (Moore, Gust et al 1984). With a more complex system, i.e. D-PS-PS-A-A, 
it has been possible to reach a lifetime of 340 psec (Gust, Moore, 1990), which means that the 


back reaction is strongly slowed. Thus one can hope to be able to use the stored chemical 
energy in later steps. 


Figure 3. The carotene-porphyrin-quinone triad 


But for electronic applications, it is necessary to incorporate these molecules in organized 
structures So that a macroscopic current could be detected. Photocurrents have been observed 
with such molecules embedded in lipid bilayers (Seta et al, 1985). Another interesting system 
has been described with A = viologen, D = ferrocene, PS = pyrene (Fujihira et al, 1985). A 


photocurrent was also observed with this molecule deposited as a L B monolayer on a 
semitransparent electrode and in contact with an electrolyte. 


14. Switch, 

If we use as bridging ligand a molecule with some "structural mobility", it is conceivable to 
control the ligand geometry by some external perturbation, and thus to control the electron 
transfer process. The simplest idea is to use a photochemical reaction. Thus several donor- 


acceptor molecules exhibit the so called TICT effect (Twisted Internal Charge Transfer) in 
which upon excitation there is a twisting motion around an essentially single bond with 


decoupling of the m systems (Grabowski et al, 1979). Thus the electronic interaction is "shut 
off" but only during the excited state lifetime. A bridging TICT molecule has been prepared 
(Launay et al, 1989). However, it should certainly be more useful to have a bistable switch by 
combining electron transfer ability and photoisomerization. But the synthetic work appears 
much more difficult. 


Another possibility to build a molecular switch is to use a complexation / decomplexation 
reaction to change the conformation of a flexible ligand. Such a system has been recently 
obtained and is presently under study in the author's laboratory (Gourdon et al, 1990). Finally 
some protonation / deprotonation reactions could be used to alter the electronic structure of the 
bridging ligand (Launay et al, 1990). 


L5. Bistable element. 


Bistable elements are already avaible if we consider photochromic molecules. Efficient systems 
have been described such as the "fulgides" (Heller, 1983) and they allow repeated cycling 
between two forms with different absorption spectra. The absorption spectra having very few 
overlap, it is possible to convert quantitatively one form into the other and to come back with a 
different wavelength (Figure 4). A special case of photochromic systems is represented by 
transition metal complexes exhibiting the spin crossover phenomenon. In some cuse they can 
also be switched between two states with different optical and magnetic properties (Hauser, 
1986). This constitutes the LIESST effect (Light Induced Excited Spin State Trapping) which is 
formally analogous to photochromism with the added dimension of magnetic effects. 


Figure 4. A bistable molecule of the"fulgide" family (Heller, 1983) 


L6. Others. 


Many other molecular devices have been proposed but have not been synthesized. One can cite 
soliton devices (Carter, F. L., 1987), a memory shift register (Hopfield et al, 1988), a logical 
gate (Aviram, 1988) A special mention should be made about molecules with easily accessible 
excited states, because they constitute potential systems to realize logical functions (see below). 


2. SOME BASIC PROBLEMS IN MOLECULAR ELECTRONICS 


Although the molecular device concept has triggered a noticeable research effort, it presently 
does not imply a revolutionary kind of experiments. Rather the synthesis and present 
characterization of these devices is made by standard chemical and physico-chemical methods. 
But for really using these molecules, it is necessary to go beyond and consider several 
problems. 


LL Testi iati ts? 


The will to devise truly unimolecular devices has led to sophisticated methods to perform 
molecular addressing, i.¢ to be able to interrogate a single molecule. This technological “tour de 
force” is presently attacked by several means. Thus an optical method to detect a single 
molecule has been described recently by Moerner (Moerner et al, 1989) and also by Orrit 
(Bernard and Orrit, 1990). Electrical addressing is already realized in STM experiments when 
imaging single molecules is achieved (Quate et al, 1989). Another electrical approach, not yet 
realized, is based on a combination of advanced nanolithographic techniques and chemical 


grafting. 


Although these ingenious methods will help the characterization of Molecular Devices and will 
pave the way for Intramolecular Physics, they cannot be considered as the ultimate goal in 
Molecular Electronics. The true challenge for building an information processing system is in 
fact to_associate a large number of components and to predict the behaviour of a large 
population. This aspect has been ignored until now because it raises formidable difficulties. The 
main one comes from the lack of a measurable intermediate quantity which would be the output 
of a device and at the same time the input of another one. In traditional electronics, this quantity 
is either the current or the potential. In Molecular Electronics the notion of electrical potential is 
not immediately apparent and one has to rely on concepts such as “quasi chemical potential” 
(Joachim, 1990, Buttiker, 1989). Another problem is to achieve the directionality in 
communication, i.c that the output of device | is the input of device 2 and not the reverse. 
Finally, since the very process of measurement introduces a quantum perturbation, there is no 
simple step-by-step procedure to build a molecular circuit. One has to concieve the overall 
system with macroscopic input and output, and give up the idea of measuring intermediate 
quantities without completely modifying the system's behaviour. 


: are 2 ' A 
Most work devoted to electron transfer has been performed with low driving forces, i.e near to 
equiliorium. This is particulary true for studies on mixed valence complexes (Woitellier et al, 
1989), and for the case of model switches embedded in a conducting chain (Sautet and 
Joachim, 1988). As a result there is!few (if any) energy dissipation in the device itself, but 
rather in the interfacing sites or the vlectron reservoirs. 


Theoretically this is not rhedibitory for an information processing system, because it has been 
shown that a computation can be performed in principle without the dissipation of energy 
(Bennett and Landauer, 1985). But it appears very difficult to put this general result in practice 
and to build an actual computing system (Keyes, 1988). As all information processing systems 
known to date proceed with energy dissipation, it is clear that we lack fundamental studies on 
the behaviour of molecular devices in the high driving force regime. Such systems would 
constitute simple 2-terminal devices, but their study is not a simple task because, if we consider 
electrically addressed molecules, a high driving force would mean an extremely high electric 
field, several orders of magnitude greater than usual values. 


As far as 3-terminal devices are concemed, the problem is to control a large flow of energy by a 
small one in order to get amplification. This aspect has not been considered until now, because 
in most cases the energy source is not clearly defined in the prototypes of molecular devices. In 
particular, it is likely that the prototypes of switches evoked in 1.4 do not realize amplification . 


23 WI Id. Molecular EI ’ in inf ‘ ; 
systems ? 


Once elaborate molecules playing the role of Molecular Devices have been obtained, it is 
tempting to imagine how to assemble them in order to realize complex functions. This line of 
thought parallels the classical development of silicium electronics from simple components 
(transistors) to more complex ones (adders, logical gates, ... ), then microprocessors and so 
on. The fact that trying to mimic closely silicium would be foolish has been recognized very 
soon, thus alternate computer architectures more adapted to molecular electronics have been 
evoked, such as cellular automata, or neural networks. However these are not miracle solutions 
because present and foreseed molecular devices can perform only very simple and rudimentary 
functions. On the other hand, they can be assembled in really great numbers, of the order of 
Avogadro number. 


To help clarify ideas, the diagram of Figure 5 can be used. Each information treatment system 
can be broadly characterized by three qualitative parameters : (i) the degree of complexity of a 
processing unit or more properly of the task it can perform (ii) the number of such processing 
units and (ili) the complexity in the connectivity between units. Thus a sequential computer (one 
very complex unit) can be represented by point SQ (see Figure 5). In parallel machines, one 
associates several processing units, but each one is still a microprocessor, only slightly 
simplified. Neural networks are presently realized with a limited number of neurons (typically 
10 -100). Their performances come from a subtle interplay between processors and synapses, 
so that the connectivity is rather complex and plays an essential réle. Finally increasing 
considerably the number of processors and simplifying both their task and the connectivity 
leads to the cellular automaton concept. Here all connexions are identical and local. 


Figure 5. Diagram showing the place of the different information processing systems 
SQ : Sequential computer (1 processor) ; NN : Neural network ; CA : Cellular automata ; 
ME : Molecular Electronics 


In this diagram, Molecular Electronics (as we can imagine it to day) would lie in the area ME 
beyond cellular automata. The complexity of the task which can be performed by a molecule is 
yery low, much lower than the one required for a cellular automaton for instance. Also the 
connectivity is very simple since chemists can only master near neighbour interactions. Only the 
number can be very great. 


Thus a computer architecture based on Molecular Electronics cannot be easily imagined at the 
present time. It does not fit in existing architectures because they have not been devised for this 
kind of components. With respect to cellular automata, the present molecules are too simple; 
with respect to neural networks, the connectivities are too rudimentary. There is clearly a need 
to “start from the bottom” , i. e. to begin to think what can be done from the point of view of 
architecture with present molecular devices. 


3. A PROGRAM FOR THE FUTURE 
ILT Lf r Hehe 


Nonlinearity appears to be 2 necessary crossing point for an information processing system. 
With a sigmoidal shape such as Figure 6 (a), one can perform threshold or logical functions. A 
shape such as shown on Figure 6 (b) yields amplification or oscillation. Both shapes can 
provide bistability if a suitable amount of feedback is introduced. Table I summarizes some 
nonlinear behaviours which can be encountered either with molecules or at the material level. 
For electrical properties, very few is known as only the behaviour at weak polarization has been 
approached (Joachim, 1990). For high polarizations one could expect resonant tunneling as in 
quantum well structures, but this is still a speculation. It is clear that experimental data on 
systems far from thermal equilibrium would be extremely desirable as they could lead to pattern 
formation (Haken, 1985). With optical properties, non linear effects are now well documented 
(Chemla and Zyss, 1987) As far as chemical effects are concerned, nonlinearities can be 
observed for the dependences of concentration with time. This could be the basis of a chemical 
parallel compute: (Haken, 1987). A demonstrator showing rudimentary image processing by a 


Figure 6. Examples of nonlinear behaviours 


chemical system has been built (Kuhnert, 1986). All these studies should be developed in a 
more systematic way. 


3.2, Logical functi ' feat aeat 


It is possible to realize logical functions with molecules. The simplest way for that is to make 
use of the properties of photochemica! excited states. Experiments in which the absorption 
spectrum of an excited state is recorded could be the basis for a logical function (AND because 
they realize a double excitation). Recent experiments on photon gated hole burning can be 
considered as the use of this AND function (Carter, T. P., Moerner et al, 1987) (see below). 
More complex functions are of course highly welcome but the realization of say a NAND gate 
would require at least three optical beams to test tie system (Birge et al, 1989). Using modified 
organic materials for non linear optics would perhaps be a realistic way to build a logical gate. It 


would be necessary to find an optical (2) or y@Q) material which could be perturbed by a 
photochemical excitation. 


Se OR ARR rE Ny . 


AO 


property single molecule material or assembly 
of molecules 


electrical strongly polarized switching materials 
molecules ’ (Potember, 1982, Denisevich, 
1981) 


optical B hyperpolarizability (2) materials 


chemical _ chemical processor 
concentration (Kuhnert, 1986) 


Table I. Summary of non linear effects. 

3.3, Cascading components 

It is vital for the future development of Molecular Electronics to study the basic association in 
which the output of a given device controls the imput of another. But for Molecular Devices, 
generally a single electron event is studied and this constitutes the signal. Thus it is necessary to 
be able to detect a single electron. With present technology it is much easier to detect a flux of 
electrons. This means that an intermediate goal would be to show how a permanent electron 
flux could be controlled by a single electron (In a later step the single electron event should 
control another single electron event). This process is already known as "Coulomb blockade" 
and occurs in ultrasmall structures when the charging energy of a single electron exceeds the 
characteristic energy of thermal fluctuations in a properly designed system (Fulton and Dolan, 
1987). The observation of this effect with a molecular device will certainly be attempted in the 
future. The use of large electron transfer molecules has been proposed (Gilmanshin and 
Lazarev, 1988) because it would theoretically facilitate the observation of single electron 
effects. 


3.4, The Ouantum Computer 


The analogy between spin glasses and neural networks leads to the idea of a quantum 
computer. This raises however considerable problems because most structures that chemists 
can build are static ones, i.e. the interactions between the different constitutive elements are 
fixed, while the essence of a neural network is to use adaptative connexions (synapses). Thus it 
is vital to investigate the mutual influences of molecular devices and to learn how to modify 
them in a dynamic way. In a similar way, the theory of Quantum Cellular Automata is just 
beginning to appear. A model has been proposed where the interactions between the different 
elernents are strictly local and are represented by a simple Hamiltonian. The time evolution has 
been computed, but only for small values of t (Grossing and Zeilinger, 1988). Thus much 
remains to be discovered. 


4. SOME POSSIBLE SHORT TERM REALIZATIONS 


As seen above, the development of Molecular Electronics still need a considerable effort in 
basic research before a complete integrated information processing sytem can be built. In the 
meantime however, there are some less ambitious but more realistic goals which are worth to be 
considered. They all have in common to be compatible with known architectures. 


AA 


LL. Synapti verial 


The implementation of neural networks suffers from the lack of a simple solution to mimic 
synapses. It would be interesting to have a material exhibiting a programmable electrical 
resistance (Hopfield et al, 1986). This effect could perhaps be realized with molecular 
materials, because there are already examples of compounds for which the conductivity can be 
altered by an electrical perturbation (Potember et al, 1982, Denisevich et al, 1981). 


2. Optical readout of electrical ial 


To overcome the "memory contention" problem which occurs in shared memory architectures, 
it has been proposed to use an optical readout of electrical potentials (Kowel et al, 1987). This 
could be performed with a thin film of molecules exhibiting strong quadratic nonlinearities, 
coated on a CCD matrix. These molecules would then modulate the intensity of a laser beam 
according to the value of the charge stored on the matrix element (electro-optic effect). 


L3.N ‘als with r feat cabl ‘bili 


Molecular materials with high cubic susceptibilities can in principle be used in optical devices 
showing bistability or four wave mixing, from which a variety of applications in optical 
treatment of signals can be conceived (Chemila and Zyss, 1987). 


4.4, Optical 


High density optical memories could be obtained by the hole burning technique which uses 
spatial and frequency addressing (Carter, T. P., Moerner et al, 1987). The variant of photon 
gated hole burning mentioned above represents a further improvement because writing an 
information necessitates two beams, so that reading can be nondestructive. Few (if any) 
molecules have been specifically designed for this kind of function, so that there is a great 
potential for improvement. 
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